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ABSTRACT 


A  comprehensive  experimental  and  analytical  Investigation  has  been 
undertaken  to  determine  the  hydroelastic  response  of  elastically-mounted 
smooth  and  rough  cylinders  In  harmonically-oscillating  flow. 

The  forces  acting  on  the  cylinder  and  the  cylinder  response  have 
been  measured  together  with  the  characteristics  of  the  flow  and  analyzed 
through  the  use  of  the  Fourier  analysis,  spectral  analysis,  Duhamel's 
integral,  and  the  wake-oscillator  model.  In  addition  the  in-line  force 
has  been  analyzed  through  the  use  of  Mori  son's  equation. 

The  results  have  shown  that  hydroelastic  oscillations  occur  in  a 
narrow  range  of  the  reduced  velocity  and  give  rise  to  significant 
amplification  of  the  exciting  force. 

The  use  of  the  Fourier  and  spectral  analysis  together  with  the 
governing  equations  of  motion  was  shown  to  predict  the  hydroelastic 
response  of  the  cylinders  in  the  synchronous  region. 

The  investigation  shed  considerable  light  on  the  interaction  between 
the  harmonic  fluid  motion  and  the  hydroelastic  response  of  cylinders  and 
provided  definitive  information  for  the  prediction  of  the  force  and 
response  characteristics  of  the  type  of  oscillations  considered  herein. 
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I.  GENERAL  SURVEY 


A.  INTRODUCTION 

Experience  has  shown  that  elastic  structures  near  linear-resonance 
conditions  can  develop  flow-induced  oscillations  by  extracting  energy 
from  the  flow  about  them.  The  oscillations  of  the  structure  modify  the 
flow  and  give  rise  to  nonlinear  interaction.  This  is  in  addition  to 
any  non-linearity  which  can  arise  from  the  restoring  force  (variable 
support  stiffness)  and/or  from  response-dependent  structural  damping. 

The  understanding  of  these  nonlinear  interactions  is  of  paramount  importance. 

Four  general  classes  cover  a  wide  range  of  the  fluid-structure 
interaction  phenomenon:  vortex- induced  oscillation,  flutter,  galloping 
and  buffeting.  This  survey  will  deal  with  only  the  first  of  these 
phenomena . 

The  problems  associated  with  vortex- induced  oscillations  and  their 
important  consequences  have  proved  to  be  exceedingly  difficult  both 
theoretically  and  experimentally.  The  difficulty  is  partly  due  to  our 
incomplete  knowledge  of  the  flow  field  about  a  stationary  body  but  it  is 
also  due  to  fundamental  problems  associated  with  the  coupling  of  the 
oscillations  and  the  flow  field.  It  would  indeed  be  gratifying  to  note 
that  our  understanding  of  the  separated  flow  about  a  bluff  body,  particu¬ 
larly  in  the  region  enclosing  the  body  and  its  near  wake,  is  complete^ 
through  laboratory  and  numerical  experiments.  This  is,  regrettably,  not 
yet  the  case,  but  there  have  been  nevertheless  significant  advances  since 
the  time  attention  has  shifted  from  idealized  vortex  streets  to  the 
details  of  separation,  three-dimensional  flow,  and  the  role  played  by  the 
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afterbody  on  the  formation,  growth,  and  motion  of  vortices.  Many  of  the 
key  developments  of  the  past  decade  have  been  incorporated  into  the 
experimental  investigation  of  flow-induced  oscillations  and  into  the 
evolution  of  conceptually  plausible  idealized  models. 

In  the  material  to  follow  emphasis  is  placed  on  the  understanding  of 
the  transverse  oscillations  of  cylinders.  No  attempt  has  been  made  to 
be  encyclopedic  in  coverage  as  there  are  recent  texts  and  survey  articles, 
for  example,  Blevins  [1],  Parkinson  [2],  King  [31,  and  McCroskey  [4], 
which  cover  the  field  in  some  detail. 

B.  VORTEX  SHEDDING  FROM  A  STATIONARY  BLUFF  BODY 

There  are  substantial  gaps  in  the  understanding  of  the  fundamental 
mechanisms  at  work  in  the  forcing  of  the  body  by  the  forces  acting  on 
its  afterbody  and  in  the  prediction  of  the  lock-in  range  in  which  resonant 
oscillations  occur.  To  a  degree,  the  difficulties  in  obtaining  a  solution 
or  in  devising  a  model  stem  from  an  inadequate  knowledge  of  the  precise 
mechanism  of  flow  separation  and  its  consequences  in  steady  flow  past  a 
stationary  body  (see,  e.g.,  Morkovin  [5],  Marris  [6],  Mair  and  Maull  [7], 
and  Berger  and  Wille  [8]  for  comprehensive  reviews  and  references). 

Observations  as  well  as  numerical  experiments  have  shown  that  the 
wake  of  a  bluff  body  is  comprised  of  an  alternating  vortex  street.  This 
phenomenon  is  intrinsic  to  the  flow  itself  and  is  a  consequence  of  the 
interaction  between  the  shear  layers,  base  pressure,  diffusion  and  dissipa¬ 
tion  of  vorticity,  and  far  wake,  i.e.,  a  consequence  of  the  particular 
vorticity  distribution  throughout  the  field. 

As  the  vortices  grow  symmetrically,  say,  in  an  impulsively-started  flow 
about  a  cylinder,  the  shear  layer  joining  the  separation  point  to  one  of 
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the  vortices  begins  to  develop  instabilities  (presumably  Tollmien- 
Schlichting  instabilities  [9])  and  is  drawn  across  the  wake  in  response 
to  the  base  pressure  reduced  by  the  action  of  the  vortex  growing  across  the 
wake  (see  Fig.  1)  [10].  This  nearly  corresponds  to  a  time  when  the  sheet 
drawn  in  has  least  circulation  or  is  most  permeable.  The  stretching, 
diffusion,  and  dissipation  of  vorticity  break  up  the  deforming  turbulent 
sheet  and  thereby  the  further  supply  of  circulation  to  the  vortex  whose 
rate  of  growth  has  already  been  reduced  to  its  minimum.  This  corresponds 
to  the  shedding  of  the  first  vortex.  The  vortex  across  the  wake  still 
continues  to  grow  (but  at  a  decreasing  rate)  and  entrains  part  of  the 
oppositely-signed  vorticity  left  in  the  wake  by  the  cut  sheet  and  the 
irrotational  fluid  drawn  from  outside  through  the  opening  created  by  the 
shedding  of  the  first  vortex.  The  shedding  process  for  the  second  vortex 
does  not  commence  until  the  circulation  in  its  feeding  sheet  decreases  to 
its  minimum,  making  the  sheet  most  susceptible  to  rapid  diffusion.  Simulta¬ 
neously,  the  sheet  deforms,  diffuses,  and  is  drawn  across  the  wake  by  the 
action  of  the  base  pressure  and  the  vortex  growing  on  the  other  side  of  the 
wake.  Then  the  shedding  cycle  repeats  itself.  The  above  mechanism, 
exhibited  by  numerical  experiments  [10]  is  quite  similar,  but  not  identical, 
to  that  suggested  by  Gerrard  [11]. 

The  separation  point  at  the  side  of  the  cut  sheet  just  passes  through 
its  minimum  angle.  The  sheet  at  the  opposite  side  of  the  wake  bears  nearly 
the  maximum  circulation,  relative  to  any  other  time,  and  its  separation-point 


For  a  circular  cylinder  there  is  a  60  degree  phase  difference  between 
the  lift  maximum  and  the  sheet  cutting  time,  i.e.,  a  vortex  is  shed  when 
the  lift  force  decreases  to  50  percent  of  the  absolute  value  of  its  maximum. 
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Initiation  of  vortex  shedding  [10] 


angle  just  passes  through  a  maximum.  The  motion  of  the  stagnation  point 
is  such  that  it  is  180  degrees  out  of  phase  with  the  separation  points, 
i.e.,  the  separation  point  is  below  the  x-axis  (flow  from  left  to  right) 
when  a  vortex  is  shed  from  the  bottom  of  the  cylinder. 

The  shedding  of  the  vortices  gives  rise  to  an  oscillating  side  thrust, 
upon  a  cylinder  of  suitable  afterbody,  in  a  direction  away  from  the  last 
detached  vortex.  The  relationship  between  the  vortex-shedding  frequency, 
cylinder  diameter,  and  the  velocity  of  the  ambient  flow  has  been  discovered 
by  Strouhal  [12]  almost  exactly  100  years  ago  in  connection  with  his  work 
on  a  special  method  of  creation  of  sound.  His  data  showed  that  the  product 
fD/V  is  nearly  constant  for  a  wide  range  of  values  of  D  and  V.  This  pos¬ 
tulated  constant  is  often  called  the  Strouhal  number  and  Strouhal 's  own 
data  suggested  a  value  of  about  0.185.  The  characterization  of  the  vortex 
shedding  process  by  a  simple  frequency  is  a  practical  simplification.  As 
first  noted  by  Lord  Rayleigh  [13],  the  Strouhal  number  is  a  function  of  the 
Reynolds  number  for  a  given  body  and  f  should  be  written  as  f  =  SQ(Re)V/D, 
(see  Fig.  2). 

At  relatively  low  subcritical  Reynolds  numbers,  the  energy  containing 
frequencies  are  confined  to  a  narrow  band,  and  the  Strouhal  number  is  about 
0.20  for  smooth  cylinders  [8].  For  Reynolds  numbers  larger  than  about 
20,000,  only  an  average  Strouhal  number  may  be  defined.  In  the  range 
0.5x10®  <  Re  <  2.5x10®  no  periodicity  is  observed  in  either  the  unsteady 


2S  =  0.21(1  -  20/Re)  for  40  <  Re  <  1,000  (Goldstein  [14]); 

S  *  0.212(1  -  12.7/Re)  in  the  range  400  <  Re  <  100,000  according  to  a 
correlation  obtained  by  Roshko  [15].  In  the  range  400  <  Re  <  1,200  the 
results  of  Roshko's  correlation  are  slightly  below  the  best-fit  line  to 
the  existing  data  [8], 
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lift  or  the  unsteady  pressures.  In  the  range  2.5x10®  <  Re  <  6.5x10® 

a  small  peak  is  observed  at  SQ  -  0.2  in  the  unsteady  lift  and  in  unsteady 

pressures  only  at  the  front  part  of  the  cylinder.  No  peak  is  found  behind 

*  6 

the  separation  point.  For  Re  >  6.5x10  ,  a  definite  peak  is  found  in  both 
the  unsteady  lift  and  the  unsteady  pressures  at  SQ  =  0.3.  The  said 
periodicity  occurs  both  at  the  front  and  rear  part  of  the  cylinder  [16]. 

The  foregoing  categorization  is  somewhat  arbitrary.  Jones  at  al., 

[17]  classified  the  unsteady  lift  force  in  three  regimes:  wide-band 
random  for  1.1x10®  <  Re  <  3.5x10®;  narrow-band  random  for  3.5x10®  <  Re  <  6x10® 
and  quasi -periodic  for  6x10®  <  Re  <  18.7x10®. 

The  spectral  content  of  the  exciting  forces  is  particularly  important 
for  bodies  which  may  undergo  in-line  and/or  transverse  oscillations  since 
the  vortex  shedding  frequency  locks  on  to  the  frequency  of  the  transverse 
oscillations  of  the  cylinder  when  there  is  a  transverse  force  and  when  the 
vortex  shedding  frequency  is  in  the  neighborhood  of  the  natural  frequency 
of  the  cylinder.  When  the  flow  about  a  freely  oscillating  cylinder  is  in 
the  critical  regime,  the  sharp  rise  in  the  Strouhal  number  does  not  occur 
and  remains  at  a  value  nearly  equal  to  that  found  at  subcritical  Reynolds 
numbers  [18]. 

Various  attempts  have  been  made  to  devise  a  universal  Strouhal  number 
which  would  remain  constant  for  differently  shaped  two-dimensional  and 
axisymmetric  bodies  (see  Table-I).  In  Table-I  h'  represents  the  wake 
width  as  determined  from  Roshko's  notched  hodograph  theory  [19];  £v, 
the  longitudinal  spacing  of  the  vortices  as  determined  by  Kronauer's 
minimum  drag  criteria  as  reported  by  Bearman  [20];  the  wake  width  at 
a  distance  of  the  wake  establishment  region;  and  ^  =  (1  -  C 
where  Cpp  is  the  base  pressure  coefficient. 
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Table  I  Universal  Strouhal  numbers 


So 

value 

Body  shape 

Source 

foh,/V><w 

0.163 

two-dimensional 

Roshko 

[19] 

wv 

0.181 

two-dimensional 

Bearman 

[20] 

foVVKw 

0.190 

axi symmetric 

Calvert 

[21] 

foVVKw 

0.163 

two-dimensional 

Simmons 

[22] 

The  foregoing  definitions  suffer  from  the  obvious  drawback  that  they 
require  the  solution  of  the  wake  formation  problem  first,  or  the 
measurement  of  one  or  more  flow  characteristics. 

Theoretical  or  semi -empirical  predictions  of  the  wake  characteristics 
have  been  attempted  by  various  researchers.  Birkhoff  [23]  demonstrated 
that  the  longitudinal  spacing  l  is  trivially  invariant  since  the  longi¬ 
tudinal  velocity  of  vortices  is  uniformly  bounded.  He  has  shown  further 
that  in  an  inviscid  fluid  the  lateral  spacing  also  remains  constant  at  a 
value  equal  to  the  initial  spacing  of  the  shear  layers.  These  led  Birkhoff 
to  the  conclusion  that  hy/i,v  is  determined  by  its  initial  value  and  not 
by  von  Karman's  stability  criteria  of  hy/Jlv  *  0.281.  The  use  of  a  wake 
oscillator  model  ("the  wake  swings  from  side  to  side,  somewhat  like  the 
tail  of  a  swimming  fish"  [23])  and  some  experimentally-determined  values 
led  Birkhoff  to  SQ  =  0.2.  One  must  keep  in  mind  that  vorticity  in  real 
vortices  is  not  concentrated  in  points,  the  vortices  are  non-circular  and 
distort  and  rotate  as  they  move  downstream  [24],  vorticity  diffuses  and 
is  swept  across  the  wake  [25],  and  dissipates  by  turbulence  18).  Vortices 
are  subjected  to  strain  fields  imposed  by  nearby  vortices.  The  resulting 
patterns  are  ever  changing  vortex  shapes  encompassing  elliptic  and 
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pearlike  geometries.  The  complexity  of  the  interaction  between  strained 
distorting  vortices  is  further  exacerbated  by  the  addition  of  turbulence 
to  the  wake,  as  this  is  likely  to  produce  a  more  diffusive  vorticity 
distribution  and  thus  an  additional  shear  field  [24].  Thus,  the  near 
constancy  of  wake  characteristics  in  the  range  of  Reynolds  numbers  where 
a  vortex  street  might  be  observed  is  primarily  due  to  the  slow  variation 
of  £v,  hy,  and  SQ  with  respect  to  the  strength  and  deformation  of  vortices. 
Nevertheless,  theoretical  idealizations  of  the  wake  give  some  clues  about 
the  asymptotic  behavior  of  the  wake.  In  particular,  one  obtains  from 
fgAy  ■  V  -  u  ,  where  V  -  ug  is  the  velocity  relative  to  the  body,  and 
from  the  use  of  von  Karman's  stability  criteria  that 

c  =  1  r/VD _  m 

0  VD  2/£4v/D)2 

where  r  represents  the  strength  of  a  vortex.  Assuming  that  D  remains 
nearly  constant  (4  / D  -  5),  one  concludes  that  the  smaller  the  strength 
of  the  vortices,  the  larger  is  the  Strouhal  number  (smaller  drag  and 
narrower  wake)  and  vice  versa.  The  vortices  which  are  subjected  to 
smaller  dissipation  linger  a  little  while  longer  in  the  near  wake  relative 
to  those  which  are  dissipated  more,  i.e.,  strong  vortices  get  stronger. 
However,  the  dependence  of  SQ  on  r/VD  is  rather  weak.  This  explains  in 
part  the  reason  for  the  success  of  the  discrete  vortex  models  in  predicting 
the  Strouhal  number  fairly  accurately  [26,27]  in  spite  of  the  fact  that 
the  calculated  vortex  strengths  are  about  35  percent  larger  than  those 
estimated  experimentally.  However  weak,  the  dependence  of  SQ  on  r/VD 
becomes  a  primary  factor  in  the  locking-on  of  the  vortex  shedding  to  the 
natural  frequency  of  the  body.  The  vortices  which  are  stronger  continue 
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to  be  fed  by  their  shear  layers  a  longer  time  period,  thus  further  reducing 
the  Strouhal  number.  In  this  process  the  mobility  of  the  separation 
points  is  important  but  not  necessary. 

Sacksteder  [28]  pursued  a  theoretical  approach  to  determine  the  Strouhal 
number  at  large  Reynolds  numbers  by  perturbing  the  d'Alembert  flow  and 
obtained  SQ  *  0.2028.  His  analysis  does  not  consider  the  wake  formation 
and  thus  it  is  not  expected  that  flows  around  objects  that  induce  large 
wakes  could  be  treated  with  a  simple  perturbation  of  the  unseparated  flow. 

The  constancy  of  the  Strouhal  number  over  a  broad  range  of  Reynolds 
numbers  does  not  imply  that  the  base  pressure  remains  constant  and  that 
a  simple  two-dimensional  vortex  emanates  from  a  separation  line.  In 
reality,  there  is  not  only  a  phase  shift  between  various  sections  of  the 
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vortex,  separated  by  a  correlation  length,  but  also  variations  in  both 
the  intensity  and  the  frequency  of  vortex  segments. 

The  variation  of  the  base  pressure  with  Reynolds  number,  in  the  range 
where  the  Strouhal  number  remains  practically  constant,  may  be  related  to 
the  variation  of  the  mean  vorticity  flux  or  to  the  variation  of  the 
correlation  length  with  the  Reynolds  number,  turbulence,  length- to-diameter 
ratio,  and  surface  roughness  [29-33].  Table-2,  as  compiled  by  King  [3], 
gives  an  approximate  idea  about  the  typical  values  of  the  correlation 
length.  The  net  effect  of  the  spanwise  variations  of  the  vortex  tube  is 
that  the  transverse  force  (lift)  coefficient  obtained  from  a  pressure 
integration  is  not  necessarily  identical  with  that  obtained  from  the  direct 
measurements  of  the  lift  force.  Partial  spanwise  correlation  leads  to 


^The  equivalent  length  over  which  the  velocity  fluctuations  in  the 
wake  may  be  described  as  perfectly  correlated. 
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Table  2  Correlation  lengths 


Reynolds  No.  range 

Correlation  length 

Source 

40  <  Re  <  150 

150  -  200 

[34] 

150  <  Re  <  105 

20  -  30 

[34] 

104  <  Re  <  4 . 5x1 Q4 

30  -  6D 

[32] 

Re  >  105 

0.050 

[34] 

Re  ■  2x1 05 

1 . 560 

[29] 

variations  in  both  the  frequency  and  the  amplitude  of  the  lift  force,  the 
variation  of  the  latter  being  more  pronounced  than  that  of  the  former. 

The  lack  of  correlation  exists  not  only  spanwise  but  also  chordwise 
[35-38]  and  the  chordwise  correlation  is  related  to  the  spanwise  correla¬ 
tion.  Comparison  between  various  results  suggests  that  with  increasing 
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Reynolds  number  over  the  range  10  to  10  ,  the  chordwise  correlation  for 
square  and  circular  cross-section  cylinders  is  maintained  or  improved.  This 
leads  to  an  increase  in  fluctuating  lift  coefficient.  The  reasons  for 
these  variations  are  not  quite  clear.  The  end  effects  [29,  39-41],  wall 
boundary  layers,  freestream  turbulence  [42-44],  and  non-uniformity  of 
the  flow  are  mentioned  often  as  possible  reasons.  The  complexity  of  the 
three-dimensional  nature  of  the  flow  about  a  cylinder  is  clearly  demon¬ 
strated  with  measurements  by  Tournier  and  Py  [45]. 

It  would  not  be  correct  to  assume  that  the  mobility  of  the  separation 
points  is  primarily  responsible  for  the  imperfect  coherence.  Even  bodies 
such  as  90-degree  wedges,  square  cylinders,  with  fixed  separation  lines, 
do  not  exhibit  perfect  correlation.  However,  the  variation  of  the  base- 
pressure  coefficient  for  bodies  with  mobile  separation  lines  is  greater 
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than  that  for  bodies  with  fixed  separation  lines  [46].  In  fact,  for 
flow  over  symmetrical  wedges  the  base  pressure  appears  to  be  insensitive 
to  the  nose  angle  [46] . 

As  noted  earlier,  the  strength  of  the  vortices  plays  an  important 
role,  particularly  in  the  near  wake.  Laboratory  and  numerical  experiments 
have  shown  that  [7,  19,  47-50]  the  net  circulation  of  a  rolled-up  vortex 
of  the  street  is  40  to  60  percent  smaller  than  that  generated  in  the  boun¬ 
dary  layer  during  a  shedding  cycle.  Prandtl  determined  that  the  initial 
vorticity  decreases  to  about  half  where  the  first  vortex  centers  appear. 
Vorticity  is  ultimately  dissipated  by  viscosity  to  which  it  owes  its 
generation.  Nevertheless,  one  may  think  of  loss  of  circulation  through 
cancellation  of  oppositely-signed  vorticity.  Primarily,  there  are  three 
mechanisms  whereby  oppositely-signed  vorticity  are  brought  close  together: 
vorticity  generated  on  the  forebody  is  carried  by  the  shear  layers  near 
that  generated  on  the  afterbody;  vorticity  of  the  deformed  and  cut  sheet 
is  carried  across  the  near  wake  by  the  entrainment  of  the  irrotational 
fluid;  and,  finally,  vorticity  is  swept  across  the  entire  wake  [25].  The 
percentages  quoted  in  the  literature  [8]  for  the  total  loss  of  vorticity 
often  imply  that  the  vortices,  once  having  acquired  a  certain  circulation, 
retain  that  circulation  throughout  the  rest  of  their  motion.  The  fact 
that  circulation  decreases  continuously  with  time  or  distance  is  demonstrated 
clearly  by  the  experiments  of  Schmidt  and  Tilmann  [50]  and  Bloor  and 
Gerrard  [48].  The  amount  of  vorticity  generated  in  the  boundary  layers 
and  the  amount  dissipated  are  of  prime  importance  not  only  for  the  flow 
about  stationary  bluff  bodies  but  also  for  those  undergoing  resonant 
oscillations.  In  fact,  the  entire  bluff  body  problem  may  be  reduced  to 
the  determination  of  the  vorticity  distribution  throughout  the  flow  field. 


This  is  not  yet  possible  for  Re  larger  than  about  100.  The  determination 
of  the  vortex  strengths  is  difficult  and  sensitive  to  the  theoretical 
and  experimental  means  employed. 

It  is  evident  from  the  foregoing  and  from  a  more  detailed  perusal  of 
the  references  cited  that  the  description  of  the  near  wake  of  a  bluff 
body  is  in  a  primitive  state.  Much  of  what  is  known  about  the  consequences 
of  separation  has  come  from  laboratory  experiments.  It  has  not  yet  been 
possible  to  develop  a  numerical  model  with  which  experiments  may  be  con¬ 
ducted  to  explain  the  observed  or  inferred  relationships  between  various 
parameters  and  to  guide  and  complement  the  laboratory  experiments.  The 
principal  difficulties  are  as  follows:  (i)  Separation  points.  They  repre¬ 
sent  a  mobile  boundary  between  two  regions  of  vastly  different  scales. 

This  in  turn  leads  to  complex  physical  nonuniformities  in  relatively 
narrow  regions  which  cannot  be  handled  within  the  framework  of  the  boundary- 
layer  theory  [51].  Finite  difference  and  Marker  and  Cell  (MAC)  techniques 
require  in  such  regions  very  small  computational  times.  The  discretization 
of  the  continuous  process  of  vorticity  generation  by  line  vortices  in  the 
vicinity  of  a  mobile  or  fixed  singular  point  (discrete  vortex  model) 
strongly  affects  the  existing  non-uniformities  and  promotes  earlier 
separation.  Attempts  to  preserve  the  prevailing  flow  conditions,  say  by 
limiting  the  influence  of  the  nascent  vortices,  while  satisfying  a  rela¬ 
tively  simple  separation  criterion  lead  to  hydrodynamical  inconsistencies 
and  non-disposable  parameters,  (ii)  Reynolds  number.  Finite  difference 
schemes  for  bluff-body  flows  are  limited  to  relatively  small  Reynolds 
numbers  whether  the  scale  of  flow  is  assumed  to  be  governed  by  a  constant 
viscosity  or  by  a  constant  eddy  viscosity.  The  large  recirculation  region 
of  the  flow  is  often  comprised  of  turbulent  vortices  even  when  the  boundary 
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layer  is  laminar.  The  transition  to  turbulence  moves  upstream  in  the  ' 

3 

shear  layers  as  the  Reynolds  number  is  increased  from  about  10  to 

4  4  1 
5x10  .  At  Re  =  5x10  ,  it  reaches  the  shoulder  of  the  cylinder  [52].  } 

i 

It  does  not  move  appreciably  further  upstream  before  the  critical  Reynolds 
number  is  reached.  Thus,  the  distribution,  turbulent  diffusion,  and 
decay  of  vorticity  and  the  interaction  between  the  wake  and  the  boundary 
layers  cannot  be  subjected  to  numerical  simulation  without  recourse  to 
some  heuristic  turbulence  models  and  inspired  foresight. 

The  representation  of  the  wake  by  clouds  of  point  vortices  or  dis¬ 
cretized  spiralling  sheets  (see,  e.g.,  [10,  26,  49,  53])  is  not  immune  to 
scaling  problems.  In  fact,  not  a  particular  Reynolds  number  but  only  a 
particular  flow  regime  may  be  specified,  depending  on  the  separation  cri¬ 
teria  used,  (iii)  Three-dimensionality.  As  noted  earlier,  even  a  uniform 
flow  about  a  stationary  cylinder  exhibits  chordwise  and  spanwise  variations. 

These  three-dimensional  effects  may  play  a  major  role  in  the  stretching  of 
vortex  filaments  and  in  the  redistribution  of  vorticity  in  all  directions. 

The  numerical  models  are  not  in  a  position  to  account  for  such  complex 
effects.  One  may  hope  to  assess  the  effects  of  three  dimensionality  by 
means  of  two-dimensional  numerical  experiments. 

It  is  against  this  background  that  the  nonlinear  fluid-structure 
interaction  projects  an  order  of  magnitude  more  complex  problems.  The 
picture  emerging  from  the  numerous  experiments  and  mathematical  models 
is  that  there  are  phenomena  to  be  explained,  parameters  to  be  quantified, 
and  conjectures  to  be  justified. 

C.  LOCK-IN  OR  SYNCHRONIZATION  IN  STEADY  FLOW 

Nunerous  experiments  have  shown  that  when  the  vortex-shedding  frequency 
brackets  the  natural  frequency  of  an  elastic  or  elastically-mounted  rigid 
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cylinder  with  a  suitable  afterbody  (capable  of  giving  rise  to  a  trans¬ 
verse  force),  the  cylinder  takes  control  of  the  shedding  in  apparent 
violation  of  the  Strouhal  relationship.  Then  the  frequencies  of  vortex 
shedding  and  the  body  oscillation  collapse  into  a  single  frequency  close 
to  the  natural  frequency  of  the  body  (see  Fig.  3).  This  phenomenon  is 
known  by  various  names:  lock-in,  locking-on,  synchronization,  hydro¬ 
elastic  or  fluidelastic  oscillations,  wake  capture,  self-controlled  or 
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self-excited  oscillations,  etc. 

The  facts  which  have  emerged  from  two  decades  of  work  on  vortex- 
induced  oscillations  may  be  summarized  as  follows: 

1.  When  a  body  is  close  to  its  linear  resonance  conditions  it  can 
undergo  sustained  oscillations  at  a  frequency  close  to  its  natural 
frequency  [55].  That  a  cylinder  should  be  excited  at  its  natural  fre¬ 
quency  when  the  frequency  of  the  exciting  force  is  equal  to  its  natural 
frequency  is  not  surprising.  But  that  the  phenomenon  encompasses  a  range 
of  ±25  to  30  percent  of  the  natural  vortex-shedding  frequency  and  that  the 
vibration  and  vortex  shedding  frequencies  lock  together  and  control  the 
shedding  process  are  surprising. 
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A  vortex-excited  oscillation  is  actually  a  forced  one  having  a 
self-excited  character  also  to  some  degree  due  to  lift  force  amplifica¬ 
tion  through  nonlinear  interactions.  By  definition,  self-excited 
oscillation  is  one  where  the  alternating  force  that  sustains  the  motion 
is  created  or  controlled  by  the  motion  itself;  when  the  motion  stops, 
the  alternating  force  disappears.  For  the  phenomenon  discussed  herein 
the  alternating  lift  force  does  not  disappear  when  there  is  no  oscilla¬ 
tion.  In  fact,  when  there  is  no  alternating  force  (a  flat  plate  normal 
to  the  flow)  there  is  no  lock-in.  Apparently,  the  often-used  definition 
of  'self-excited  oscillation'  is  a  misnomer. 
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Fig.  3  Response  characteristics  of  a  freely-oscillating 
cylinder,  from  Feng  [54],  (see  also  [2]) 


2.  The  interaction  between  the  oscillation  of  the  body  and  the 
action  of  the  fluid  is  nonlinear  (e.g.,  [56-57]). 

3.  The  vortex  shedding  does  not  necessarily  result  in  an  alternating 
transverse  force.  There  must  be  a  suitable  afterbody,  and  hence  an 
alternating  lift  force,  and  the  shedding  frequency  should  coincide  with 
or  be  near  the  natural  frequency  of  the  body,  or  one  of  its  harmonics, 
for  sustained  oscillations  to  occur.  Thus,  the  magnitude  and  occurrence 
of  sustained  oscillations  strongly  depend  on  the  lift  coefficient  of  the 
stationary  body. 

4.  The  parameters  Vr  =  V/ fnD  and  K.  =  mc/pD2  [58]  or  Rp=  mc/pLD^C®)^ 
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[60]  or  Ar  =  (27t^) (ttSo)  /pr  are  of  major  importance  in  determining  the 
amplitude  of  oscillations  and  the  range  of  synchronization  for  a  given 
body  [59-61]. 

5.  The  velocity  range  over  which  synchronization  occurs  increases 
with  oscillation  amplitude  [62]. 

6.  The  excitation  range  of  cross-flow  oscillations  in  air  extends  over 
4.75  <  Vr  <  8  and  maximum  amplitudes  occur  in  the  range  5.5  <  Vr  <  6.5 

[3].  In  water,  the  excitation  range  of  the  transverse  oscillations  can 
be  increased  to  4.5  <  Vp  <  10  with  maximum  amplitude  falling  within  the 
range  of  6.5  <  Vr  <  8. 

7.  For  a  circular  cylinder  with  large  L/D,  synchronization  begins 

when  f  n  f  and  ends  at  about  f./f  pp  1.4.  The  maximum  amplitude  occurs 
on  on  r 

near  the  middle  of  this  range.  At  the  end  of  the  lock-in  range,  vortex- 
shedding  frequency  jumps  to  that  governed  by  the  Strouhal  relationship, 
but  the  cylinder  continues  to  oscillate  at  f  2.  f  .  This  is  true  at 
both  ends  of  the  lock-in  range  [54]  and  shows  that  the  response  is  not  a 
simple  forced  vibration  at  the  exciting  natural  Strouhal  frequency. 
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8.  For  a  cylinder  with  fixed  separation  points  (e.g.,  a  D-shaped 

cylinder)  synchronization  extends  over  a  range  of  0.96  to  1.1  of  the 

resonant  velocity  (f  /fp  from  about  0.85  to  1.1)  [24].  According  to 

Feng  [54],  the  said  range  extends  over  a  range  of  5  <  Vr  <  7.  Gowda  [39] 

found  that  the  lock-in  begins  already  when  fQ  is  about  80  percent  of  f 

and  ends  when  f  r.  f„.  The  maximum  amplitude  occurs  near  the  end  of  the 
on 

synchronization. 

9.  In  synchronization,  the  slantwise  vortex  shedding  is  replaced  by 
parallel  vortex  shedding  [44,  62]. 

10.  The  correlation  length  increases  rapidly  with  amplitude  [57,  62-65], 
The  increase  of  the  correlation  length  in  smooth  flow  is  much  larger  than 

in  turbulent  flow  [66].  In  smooth  flow  the  correlation  length  is  estimated 

by  numerical  extrapolation  to  increase  from  about  3.5D  to  400  for 
4  5 

Re  =  1.9x10  in  the  range  0.05  <  A/D  <  0.1.  In  turbulent  flow,  it  is 
again  estimated  to  vary  from  about  2. 50  to  10D  in  the  same  A/D  range. 

The  rate  of  increase  is  steeper  than  linear  but  does  not  show  any  abrupt 
change  which  would  indicate  a  sudden  development  of  the  lock-in  once  a 
threshold  amplitude  is  achieved. 

11.  For  bodies,  with  mobile  or  fixed  separation  points,  undergoing 
sustained  oscillations,  the  vortex  strength  is  increased  [24,  67].  This 
could  be  either  or  both  due  to  increased  rate  of  vorticity  flux  or  due 
to  less  destruction  of  vorticity  in  the  near  wake.  The  growing  vortex 
on  a  vibrating  cylinder  seems  to  roll  up  more  quickly  [24]. 

5 

The  generality  of  these  values  is  not  certain  since  in  Novak  and 
Tanaka's  [66]  experiments  the  cylinder  length  (between  the  end  plates) 
was  only  12  diameters. 


12.  The  transverse  force  needed  to  excite  a  cylinder  to  large  ampli¬ 
tude  oscillations  is  far  greater  than  that  exerted  by  vortex  shedding, 
i.e.,  the  transverse  force  is  amplified  (e.g.,  [54],  [56],  [61]).  The 
increase  of  the  correlation  length  is  responsible  only  for  about  10 
percent  of  the  lift  increase,  i.e.,  it  is  not  the  major  factor  for  the 
lock-in  phenomenon. 

13.  A  hysteresis  behavior  may  exist  in  the  amplitude  variation  and 
frequency  capture  depending  on  the  approach  to  the  resonance  range  -- 
whether  from  a  low  velocity  or  from  a  high  velocity  [54].  No  universal 
behavior  is  noted  and  the  reasons  are  not  clear.  The  jump  condition 
(double  amplitude  response)  may  originate  in  the  fluid  system,  and  there¬ 
fore  in  the  lift  force,  and  not  in  the  cylinder  elastic  system  [2],  The 
'combination-oscillation'  model  [68],  (based  on  the  assunption  that  the 
lift  oscillator  has  two  components,  one  at  the  cylinder  frequency,  repre¬ 
senting  lock-in,  and  the  other  approximately  at  the  frequency  of  vortex 
shedding  for  the  stationary  cylinder)  gives  rise  to  jumps  between  two 
stable  oscillation  amplitudes  but  only  at  the  ends  of  the  lock-in  range 

and  not  in  the  middle  of  the  said  range,  as  found  by  Feng  [54],  (see  Fig.  3). 
The  jump  may  be  the  consequence  of  a  variable  structural  damping  (dependent 
on  the  interaction  of  the  ambient  velocity  and  structure  support,  vibra¬ 
tion  velocity  and/or  amplitude,  etc.)  [69-70]  or  of  a  nonlinear  spring 
behavior  [71].  It  is  also  entirely  possible  that  the  jump  condition 
originates  in  both  the  fluid  and  the  cylinder  elastic  system,  including 
blockage  effects  [72]  and  the  vibrations  of  the  test  apparatus  [54] . 

14.  The  absolute  value  of  the  base  pressure  and  the  in-line  or  drag 
force  increase  with  the  amplitude  of  oscillation  [3,  56,  61,  72,  73].  For 
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a  circular  cylinder,  Cp^  increases  as  much  as  60  percent  for  A/D  =  0.29 
[72].®  For  a  D-shaped  cylinder,  Davies  [24]  found  an  increase  of  70 
percent  relative  to  the  stationary  D-cylinder  value. 

15.  The  phase  difference  between  the  displacement  and  force  signals 
increases  almost  linearly  at  first  and  then  rapidly  [54,  56]. 

16.  The  mobility  of  the  separation  points  is  important  but  not 
necessary  for  synchronization.  The  separation  point  on  a  circular  cylinder 
moves  over  an  arc  whose  magnitude  depends  on  the  frequency  and  amplitude 
of  the  oscillation.  The  magnitude  of  the  separation-point  travel  reaches 

a  peak  when  the  cylinder  frequency  is  about  0.9  of  the  wake  frequency 
[74,  75]. 

17.  The  actual  instantaneous  value  of  the  wake  angle  is  greater  than 
that  between  the  ambient  flow  velocity  and  the  relative  fluid  velocity 
[74,  75]. 

18.  The  overall  wake  width  is  practically  unaffected  in  spite  of  the 
lateral  movement  of  the  body  that  generates  it.  Koopman  [62]  inferred 
from  flow-visualization  studies  an  apparent  decrease  in  lateral  vortex 
spacing  in  the  resonant  wake  of  an  oscillating  circular  cylinder.  Davies 
[24]  found  about  a  20  percent  increase  for  a  D-shaped  cylinder  synchronized 
in  a  flow  at  much  high  Reynolds  numbers  (7x10  <  Re  <  4x10  ,  A/D  =  0.2). 
Griffin  [67]  found  that  the  lateral  spacing  of  the  vortices  decreases 

with  increasing  amplitude  and  is  unaffected  by  changes  in  frequency  of 
oscillation. 

^Fixing  separation  at  0  a  t 70  degrees  by  a  tripping  wire  on  a  circular 
cylinder  yielded  a  nearly  constant  C  .  at  -1.39  for  A/D  =  0.1  for  all 
values  of  f/f  from  0.4  to  1.3  [72]?°  The  reason  for  this  anomalous 
behavior  1scnot  clear.  Confirmatory  experiments  are  needed. 


19.  Longitudinal  spacing  does  not  change  with  the  amplitude  of 
motion  in  the  synchronization  region  but  it  varies  inversely  with  the 
frequency  ratio  [77].  For  a  cylinder  undergoing  forced  transverse 
oscillations  with  the  frequency  f  ,  £v  decreases  from  its  value  at 

f  *  fR  as  f  /fQ  increases  and  vice  versa.  In  other  words,  the  wake 
contracts  for  f  /f„  >  1  and  expands  for  f /f  <  1 . 

20.  The  formation  length  l ^  exhibits  a  similar  variation  with  f  /fQ. 
However,  unlike  the  longitudinal  spacing,  the  formation  length  decreases 
systematically  with  increasing  amplitude  of  vibration  [76].  For  a  given 
f  /f  ,  the  vortices  in  the  near  wake  go  through  a  spacing  adjustment  from 
an  amplitude-dependent  spacing  to  an  ampl itude- independent  longitudinal 
spacing. 

21.  Forced  vibrations  may  differ  from  the  naturally  occurring  vortex- 
induced  vibrations  since  the  latter  depend,  to  some  extent,  on  the  history 
of  the  fluid  motion.  History  effects  may  be  particularly  important  in 
cases  where  hysteresis  occurs.  Oriven  models  overshadow  the  intricate 
interaction  between  the  evolution  of  the  motion  and  the  variations  of 

the  pressure  distribution. 

22.  Frequency  demultiplication  or  frequency  division  (existence  of 
the  natural  vortex-shedding  frequency  in  the  state  of  synchronization 
in  the  wake)  may  occur  when  the  forced  cylinder  frequency  f  is  an 

V» 

integral  multiple  of  the  shedding  frequency^  [56,  57,  72,  77]. 

7  -k 

For  example  x  =  2b  ^sin  t  is  the  exact  solution  of 

x  -  2au(l-bx^)x  +  cd^x  *  -4ab‘Jssin(3ajt) 

and  is  a  special  case  of  the  phenomenon  called  'frequency  demultiplication' 
by  van  der  Pol  [783 . 
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23.  In-line  oscillations  occur  within  two  adjacent  regions  [79-86]. 

The  first  is  in  the  range  of  1.25  <  Vr  <  2.5,  maximum  amplitudes  occurring 

at  V„  =  2.1.  The  second  region  extends  from  V„  =  2.7  to  1/  =  3.8  with 

r  r  r 

maximtm  amplitudes  at  Vr  s  3.2.  The  first  instability  region  is  accom¬ 
panied  by  symmetric  vortex  shedding  (as  if  the  flow  started  impulsively 
from  rest  at  each  cycle)  and  the  second  by  alternate  vortex  shedding. 

Several  mathematical  models  have  been  proposed  [87-91],  (see  also 
[2,  55,  71,  and  92]  for  additional  models)  in  an  attempt  to  simulate 
and/or  explain  some  of  the  experimental  results  summarized  above.  These 
models  do  not  include  the  analysis  of  the  flow  field  and  the  fluid- 
mechanical  arguments  invoked  in  their  evolution  are  not  altogether  con¬ 
vincing.  Thus,  their  worth  should  be  measured  not  so  much  by  their 
capacity  to  obtain  functional  relations  among  significant  parameters 
that  lead  to  the  basic  understanding  of  the  phenomenon  but  rather  by 
their  ability  to  produce  results  which  are  qualitatively  similar  to  those 
obtained  experimentally. 

The  most  noteworthy  among  the  oscillator  models  is  the  one  proposed 
by  Hartlen  and  Currie  [87]  where  a  van  der  Pol  type  soft  nonlinear 
oscillator  for  the  lift  force  is  coupled  to  the  body  motion  by  a  linear 
dependence  on  cylinder  velocity.  This  model  is  based  partly  on  a  sugges¬ 
tion  by  Birkhoff  and  Zarantonello  [23],  mentioned  earlier  in  connection 
with  the  estimation  of  the  Strouhal  number,  and  by  Bishop  and  Hassan  [56] 
in  connection  with  their  now  classic  experiments  with  oscillating  cylinders 
in  uniform  flow.  The  model  has  its  roots  in  mechanics  and  electricity 
rather  than  in  the  equations  of  fluid  motion.  The  details  of  this  model 
will  be  described  later. 
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0.  LOCK  IN  OR  SYNCHRONIZATION  IN  OSCILLATING  FLOW 

There  has  been  relatively  little  theoretical  and  experimental  work 
on  the  complex  dynamic  response  of  elastic  structures  to  oscillating 
flows. 

Laird  [93]  explored,  in  1962,  the  effects  of  support  flexibility  by 
oscillating  a  vertical  cylinder  through  still  water.  He  found  that  the 
forces  acting  on  a  flexibly-supported  oscillating  cylinder  can  exceed 
4.5  times  the  drag  force  on  the  cylinder  rigidly-supported  while  moving 
at  a  uniform  velocity  equal  to  the  maximum  velocity  during  the  oscillation 
and  that  a  cylinder,  flexible  enough  to  have  transverse  oscillations  with 
amplitudes  more  than  half  the  diameter,  while  performing  large  amplitude 
oscillations  in  water,  tends  to  oscillate  transversely  at  the  eddy-shedding 
frequency  and  to  vibrate  at  twice  the  eddy  frequency  in  the  in-line 
direction.  Laird's  work  pointed  out  the  significance  of  the  vibrations 
in  harmonic  flow  and  the  need  for  a  detailed  investigation  of  the  phenomenon. 

Vaicaitis  [94]  investigated  the  response  of  deep-water  piles  due  to 
cross-flow  forces  generated  by  wind-induced  ocean  waves.  The  resulting 
cross-flow  forces  were  treated  as  random  processes  in  the  time-space 
domain  and  are  assumed  to  be  dependent  on  fluid  velocities  and  vortex 
shedding  processes.  Out  of  necessity>  Vaicaitis  had  to  make  a  number  of 
assumptions  regarding  the  structural  motion,  the  added  mass  coefficient, 
Strouhal  number,  etc.  Vaicaitis  was  careful  to  note  that  some  of  his 
conclusions  were  based  on  a  particular  example  and  that  they  might  not 
be  true  for  piles  with  different  dynamic  characteristics. 

Selna  and  Cho  [95]  investigated  the  in-line  resonant  response  of  a 
tall  structure  by  assuming  constant  drag  and  inertia  coefficients  in 
calculating  the  exciting  fluid  force  through  Morison's  equation  [96], 
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modified  to  take  into  account  the  motion  of  the  structure.  Their  calcu¬ 
lations  have  shown  that  a  resonant  motion  is  possible  and  that  the  dynamic 
deflection  of  a  flexible  structure  can  considerably  exceed  its  static 
deflection. 

Verly  and  Every  [97]  measured  wave-induced  stress  on  similar  rigid 
and  flexible  vertical  cylinders  in  a  wave  channel  at  relatively  low 
Keulegan-Carpenter  and  Reynolds  numbers.  Even  though  they  were  unable 
to  correlate  their  data  with  any  suitable  parameter  governing  the  motion, 
they  concluded  that  the  vibration  is  caused  by  the  cylinder's  response  to 
eddy  shedding  and  that  there  is  no  fluid-structure  interaction.  They 
found  that  the  vibration  occurs  if  (Um  represents  the  maximum 

velocity  of  flow  or  wave  in  a  cycle;  and  D,  the  diameter  of  the  cylinder) 
is  greater  than  about  unity  for  any  natural  frequency,  wave  frequency, 
and  damping.  The  Um/fnD  parameter  never  reached  high  enough  values  in 
Verly  and  Every' s  experiments  for  the  cylinder  to  undergo  synchronized 
oscillations  as  the  present  investigation  has  shown. 

McConnell  and  Wilson  [98]  conducted  experiments  similar  to  those  of 
Laird  [93]  by  sinusoidally  driving  a  circular  cylinder  in  still  water. 
These  experiments  were  conducted  at  rather  small  Reynolds  numbers  with 
0.25  in  and  0.50  in  cylinders.  Furthermore,  these  investigators  measured 
only  the  amplitude  of  oscillation  and  calculated  the  lift  coefficient  at 
resonance  from  a  slightly-modified  version  of  an  equation  developed  by 
Sarpkaya  [60],  They  have  concluded  that  Um/fnD  is  one  of  the  most  impor¬ 
tant  parameters  governing  the  resonant  oscillations. 

Sawaragi  et  al.,  [99]  investigated  the  in-line  and  transverse  dynamic 
response  of  a  cantilevered  circular  cylinder,  with  a  concentrated  mass 
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at  its  top,  in  waves  of  small  amplitude.  The  Reynolds  number  ranged 
from  1,500  to  6,200  and  the  rms  value  of  the  Keulegan-Carpenter  number 
(K  *  UmT/D,  calculated  over  the  submerged  length  of  the  cylinder  through 
the  use  of  rms  value  of  the  maximum  of  the  horizontal  velocities)  ranged 
from  2  to  about  20.  They  have  approximated  the  lift  coefficient  for  a 
rigid  cylinder  by  a  Rayleigh  distribution  and  calculated  the  dynamic 
response  of  the  test  cylinder.  The  results  are  of  limited  value  since 
the  interaction  between  the  synchronization  and  the  force  amplification 
was  ignored  and  since  no  more  than  three  vortices  were  shed  during  a  half 
cycle.  Furthermore,  the  rapid  variation  of  the  Reynolds  number,  Keulegan- 
Carpenter  number,  vortex-shedding  coherence,  and  lift  force  with  depth, 
in  a  range  of  Keulegan-Carpenter  numbers  where  the  phenomenon  is  least 
understood,  obscured  the  role  played  by  the  governing  parameters. 

It  is  clear  from  the  foregoing  that  there  are  no  well-established 
analytical  methods  either  to  analyze  the  data  or  to  predict  the  character¬ 
istics  of  the  flow  or  of  the  vibration.  Thus,  it  was  decided  to  undertake 
a  comprehensive  investigation  of  the  dynamic  response  of  elastically- 
mounted  smooth  and  rough  cylinders  in  harmonic  flow  particularly  at  high 
Reynolds  numbers  and  Keulegan-Carpenter  numbers. 
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II.  EXPERIMENTAL  EQUIPMENT  AND  PROCEDURES 


A.  U-SHAPED  OSCILLATING-FLOW  TUNNEL 
1 .  Original  Design 

Experiments  carried  out  by  Sarpkaya  and  Tuter  [TOO]  on  small 
smooth  cylinders  at  low  Reynolds  numbers,  in  1974,  showed  the  necessity 
of  an  apparatus  capable  of  producing  larger  Reynolds  numbers.  Various 
designs  were  considered,  ranging  in  a  broad  spectrum  of  complexity  and 
versatility.  One  of  the  options  was  to  oscillate  the  fluid  by  a  suitable 
piston  mechanism.  The  achievement  of  desired  amplitudes  and  frequencies 
needed  to  obtain  high  Reynolds  and  Keulegan-Carpenter  numbers  required  an 
extremely  complex  system.  Furthermore,  the  vibrations  generated  by  the 
piston  motion  could  not  have  been  easily  overcome.  Thus,  this  design  has 
been  abandoned.  Finally,  among  various  options,  a  U-shaped  oscillating- 
flow  tunnel  was  designed  and  built  [101].  The  tunnel  in  its  original 
design  consisted  of  11  modules  for  ease  of  construction,  transportation, 
and  mounting.  Each  module  is  made  of  3/8  inch  aluminum  plates  and  rein¬ 
forced  with  1/2x4x18  inch  aluminum  flanges  welded  to  the  plates.  The 
modules  were  assembled  with  the  help  of  an  air  dyring  silicon-rubber  seal 
between  the  flanges  of  two  adjacent  modules  and  1  inch  steel  bolts  placed 
6  inches  apart.  The  inside  of  each  module  was  precision  machined  so  that 
the  largest  misalignment  was  about  0.04  inches. 

The  cross-section  of  the  two  legs  is  6x3  feet  whereas  that  of  the 
test  section  is  3x3  feet.  This  selection  was  dictated  by  several  factors 
such  as: 

(a)  available  ceiling  height; 

(b)  pressures  to  be  encountered; 


(c)  desire  to  obtain,  at  the  test  section,  an  amplitude  or 
velocity  of  oscillation  at  least  twice  that  of  the  free 


surface; 

(d)  period  of  oscillation; 

(e)  desired  Reynolds  numbers  and  relative  amplitudes; 

(f)  possible  separation  of  the  flow  at  the  corners; 

(g)  natural  damping  of  the  oscillation;  and 

(h)  magnitude  and  the  frequency  of  the  forces. 

The  length  of  the  horizontal  test  section  was  chosen  larger  than  twice 
the  actual  amplitude  (amplitude  at  the  test  section  or  virtual  amplitude) 
to  ensure  fully  developed  uniform  flow  at  the  test  section.  Finally,  the 
two  corners  of  the  tunnel  were  carefully  streamlined  to  prevent  flow 
separation  (see  Fig.  4).  This  design  proved  to  be  more  than  adequate  for 
the  experiments  planned  at  that  time. 

The  auxiliary  components  of  the  tunnel  consisted  of  plumbing  for 
the  filling  and  emptying  of  the  tunnel  with  hot  and  cold  water  (50°  to 
120°F),  a  heat  exchanger,  butterfly-valve  system,  and  the  air  supply 
system.  The  tunnel  capacity  is  about  6000  gallons  and  could  be  filled 
by  city  water  or  by  means  of  a  pump  from  a  reservoir. 

The  butterfly-valve  system  (mounted  on  top  of  one  of  the  legs  of 
the  tunnel)  consisted  of  four  plates,  each  18  inches  wide  and  36  inches 
long.  A  1  inch  steel  shaft  was  placed  at  the  axis  of  each  valve  plate 
(see  Fig.  5).  Aluminum  housings  supported  both  ends  of  the  shaft  with 
self-aligning  ball  bearings.  A  6  inch  gear  was  attached  to  one  end  of 
each  shaft  which  extended  beyond  the  bearing.  All  four  valve  plates  were 
then  aligned  and  driven  by  a  simple  rack  and  piston  system.  The  rack 
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Schematic  drawing  of  the  initial  version  of  the  U-shaped  water  tunnel 
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as  actuated  by  an  air-driven  piston  with  the  help  of  a  three-way  valve 
connected  to  the  laboratory  air  supply  system. 

To  set  the  fluid  in  harmonic  oscillation,  butterfly  valves  at 
the  top  of  one  of  the  vertical  legs  were  closed  allowing  the  sealed 
chamber  to  be  charged  with  high  pressure  air.  This  above-ambient  pressure 
caused  the  water  tunnel  to  rise  in  the  opposite  leg,  left  open  to  the 
atmosphere.  The  three-way  control  valve  enabled  one  to  release  the  air 
quickly  by  sudden  opening  of  the  butterfly- val ve  and  put  the  bulk  of 
fluid  into  a  harmonic  oscillation.  The  damping  of  the  fluid  oscillation 
was  in  a  very  acceptable  range  in  the  case  for  the  experiments  with  the 
fixed  cylinders.  But,  as  will  be  seen  later,  design  improvement  was 
necessary  in  the  case  of  oscillating  cylinders. 

2.  First  Design  Improvement 

In  order  to  investigate  the  effect  of  cylinder  oscillations  on 
the  transverse  forces  acting  on  it,  it  was  necessary  to  construct  an 
elastic  mounting  with  only  one  degree  of  freedom.  Several  design  options 
were  considered.  Finally,  the  requirements  for  rigidity,  strength  and 
torsional  stability  in  the  in-line  direction  quickly  led  to  the  use  of 
helical  steel  springs.  In  the  final  design,  the  cylinder  was  rigidly 
fixed  to  two  low-friction  linear  bearings  resting  on  two  matched  helical 
springs.  The  motion  of  the  bearings  was  constrained  by  two  rods  perpen¬ 
dicular  to  the  direction  of  flow  as  shown  in  Fig.  6. 

The  force  transducers  were  mounted  as  shown  in  Fig.  7  to  the 
aluminum  split-bearing  housing  blocks.  Each  block  was  manufactured  to 
accommodate  an  adjustable  ball  bearing  with  a  nominal  2.5  in.  bore  and  a 
length  of  5.0  in.  The  bearings  rode  on  2.5  in.  diameter  stainless  steel. 
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hardened  and  ground,  solid  shafts^  Adji^tmertfs  of„th*  bearing  d«<metar  * 
eliminated  bearing  chatter.  The  springs  were  centered  around  the  shafts 
by  a  retaining  ring  on  the  housing  blocks  and  by  the  yoke  support  at  the 
bottom. 

The  yoke  supports  allowed  the  blocks  and  springs  to  be  raised  or 
lowered  to  center  the  cylinders  vertically  in  the  test  section  and  per¬ 
mitted  the  use  of  various  cylinders,  springs  and  weights  while  accommo¬ 
dating  the  effects  of  the  buoyant  force.  Each  yoke  was  positioned  by 
another  yoke  which  was  positioned  with  precision  by  a  threaded  third  brass 
rod  passing  through  the  upper  yoke  and  lifting  against  the  housing  top. 
Alignment  of  mounts  and  rods  was  accomplished  by  a  jig  which  was  placed 
between  the  two  blocks,  leveled,  then  raised,  and  lowered  to  position  the 
two  stainless  steel  rods.  Detailed  drawings  of  the  system  are  shown  in 
Figs.  8  and  9. 

Figure  10  shows  a  dimensioned  drawing  of  one  of  the  two  support 
housings.  Each  was  welded  of  3/8  in.  aluminum  plate.  The  windows  were 
covered  by  removable  Plexiglas  plates  for  access  for  calibration  and 
maintenance,  or,  when  mounted  in  place,  for  the  observation  of  cylinder 
movement  during  tunnel  operation.  Vents  and  drains  were  installed  at 
top  and  bottom  of  the  housing  to  remove  air  or  water  in  filling  or  draining. 
The  housings  were  flanged  and  bolted  to  the  U-tunnel  wall  with  0-ring 
seals  for  water  proofing  as  may  be  seen  in  Fig.  9. 

A  6.0  in.  circular  access  hole  was  cut  in  the  tunnel  walls  at 
the  test  section.  Each  hole  was  covered  with  a  plate  of  identical  size. 

A  slot  of  1  in.  by  6  in.  was  cut  in  each  cover.  The  covers  allowed  the 
installation  of  the  cylinders;  and  the  slots,  the  unobstructed  linear 
oscillation  of  the  force-transducer  arms  attached  to  the  cylinder. 
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Fig.  8  Side  view  of  a  cylinder-support  mount  and  housing 
(dimensions  in  cm) 


The  usefulness  of  this  design  stems  from  the  fact  that  the 
support-housing  assembly  could  be  rotated  90  degrees,  thus  enabling  one  to 
conduct  experiments  for  cylinder  oscillations  in  the  in-line  direction. 

3.  Second  Design  Improvement 

A  body  oscillating  in  a  fluid  will  extract  energy  from  it. 
Therefore  one  might  expect  this  energy  loss  would  add  up  to  the  natural 
damping  of  fluid  oscillation  in  the  tunnel  giving  rise  to  a  faster  rate 
of  decay  of  fluid  oscillation  amplitudes.  This  would  mean  that  lock-in 
phenomenon  could  not  be  observed  over  many  cycles  of  oscillation.  Experi¬ 
ments  performed  by  Klinkhamer  [102]  suggested  the  necessity  of  devising  a 
system  capable  of  restoring  the  energy  lost,  due  to  damping,  to  the  fluid. 
This  would  enable  one  to  maintain  the  amplitude  of  the  flow  oscillation 
at  a  desired  level  thus  allowing  the  observation  of  the  lock-in  over  a 
large  number  of  flow  cycles.  For  this  purpose  an  air-supply  system  was 
designed  to  introduce  large  quantities  of  air  in  one  of  the  legs  of  the 
tunnel  during  prescribed  time  intervals.  The  necessary  air  was  supplied 
by  a  centrifugal  fan  driven  by  a  3HP  electric  motor.  The  fan  was  first 
connected  to  a  30x30x33  inch  air-control  box  which  regulated  the  quantity 
of  air  supplied  to  the  tunnel.  The  air-control  box  consisted  of  a  gate 
valve  at  the  exhaust  of  the  fan  and  another  gate  valve  on  top  of  the  box. 
The  first  gate  valve  was  raised  and  lowered  vertically  by  a  threaded  rod 
attached  to  a  manually  operated  gear.  This  allowed  a  fine  adjustment  of 
the  air  supply.  The  gate  valve  on  top  of  the  box  was  activated  by  moving 
it  horizontally  by  hand  for  a  coarse  adjustment  of  the  air  supply. 

The  air  was  introduced  into  the  tunnel  by  a  duct  through  a  simple 
control  system  (see  Fig.  11).  The  function  of  the  control  system  was  to 
introduce  the  air  into  the  tunnel  at  prescribed  time  intervals  so  that 
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1.  Blower. 

2.  Air-control  box. 

3.  Vertical  gate  valve. 

4.  Horizontal  gate  valve. 

5.  Air  duct. 

6.  Air  valve. 

7.  Adjustable  connecting  rod. 

8.  Sprocket  and  chaine. 

9.  Variable  speed  motor. 

10.  Flywheel. 

11 .  Butterfly  valve. 

12.  Vertical  leg  of  the  tunnel. 


Fig.  11  Sketch  of  the  air  supply  system 


the  energy  lost  due  to  damping  could  be  restored  to  the  fluid  (see  Fig. 

12).  The  air  control  system  con^stpd  ofj  drum^af  2«£eeW'n  4’iame'ter*  **■ 
and  3  feet  long.  Three  openings  at  angles  as  specified  in  Fig.  12  were 
cut  into  the  drum.  The  door  number  (1)  functioned  as  inlet  to  the  tunnel 
and  door  number  (2)  worked  as  exhaust  to  the  ambient.  Inside  the  drum 
a  vane  of  the  same  size  as  the  diameter  and  length  of  the  drum  was  placed. 
The  vane  was  rocked  about  a  mean  position  by  an  adjustable  connecting 
rod  through  a  sprocket  wheel  and  chain  mechanism.  The  system  was  driven 
by  a  variable  speed  electric  motor  and  gear.  Two  flywheels,  one  attached 
to  the  drive  shaft  of  the  connecting  rod  and  the  other  to  the  motor  shaft 
assured  smooth  and  uniform  rotation.  An  RPM  counter  switch  was  installed 
at  the  other  end  of  the  drive  shaft  and  the  signal  was  electrically  trans¬ 
mitted  to  a  digital  counter  on  the  instrument  panel. 

The  top  of  the  leg  of  the  tunnel  to  which  air  was  supplied  had  to 

be  sealed  to  prevent  air  from  escaping.  This  was  done  by  plywood  boards 

in  such  a  way  so  as  to  create  an  air-tight  chamber  on  the  top  of  the 
tunnel  leg. 

With  this  system  very  high  amplitudes  of  flow  oscillation  were 
achieved.  This,  in  turn,  necessitated  the  lengthening  of  the  horizontal 
section  of  the  tunnel  in  order  to  maintain  uniform  flow  about  the  test 
cylinders.  For  this  purpose  two  additional  modules  (each  2.5  ft  long)  of 
the  same  cross  section  as  the  existing  test  section  were  added  to  each 
side  of  the  test  section.  On  both  sides  of  these  new  modules^  two  side 

windows  (14x19  in.)  were  cut.  One  of  the  modules  also  had  a  door  on  top 

of  it.  The  side  windows  were  covered  with  plexiglass  for  the  observation 
of  the  cylinder  motion.  However  the  main  purpose  of  these  windows  was 
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to  provide  an  easy  and  quick  access  into  the  tunnel  for  the  mounting  of 
the  cylinders  and  Tore?  calibrations.  Finally,  a  new  water  level  sight 
window  was  cut  into  the  wall  of  the  vertical  leg  to  facilitate  the  reading 
of  the  water  level.  The  old  butterfly-valve  system  on  top  of  the  other 
leg  was  not  removed  since  by  closing  it  the  flow  oscillation  could  be 
quickly  stopped. 

Figures  13a  and  13b  show  photographs  of  the  two  sides  of  the  tunnel 
after  its  completion.  Figure  14  shows  a  schematic  drawing  of  the  tunnel. 

The  natural  period  of  the  tunnel  after  these  modifications  turned 
out  to  be  5.960  seconds. 

To  put  the  fluid  in  the  tunnel  into  oscillation,  the  RPM  of  the 
variable  speed  motor  driving  the  air-valve  was  set  at  such  a  value  that 
the  period  of  the  rotation  of  the  vane  was  exactly  the  same  as  that  of 
the  tunnel.  Once  the  fan  was  turned  on,  every  time  the  vane  cleared 
port  (1)  of  the  air-valve,  air  flowed  into  the  tunnel  and  pushed  the  water 
down.  The  time  of  opening  of  the  port  (1)  was  exactly  synchronized  to 
the  time  when  the  water  had  reached  its  maximum  amplitude  and  began  to 
fall  down.  The  pressure  of  the  supplied  air  restored  the  necessary 
energy  to  the  water  to  maintain  its  amplitude  at  a  certain  level.  By 
opening  or  closing  the  gate-valves  on  the  air-control  box,  one  could 
change  the  quantity  of  air  and  consequently  the  amplitude.  This  design 
turned  out  to  be  extremely  successful  because  the  tunnel  could  sustain 
oscillations  at  a  desired  amplitude  for  an  indefinite  period  of  time. 

B.  CIRCULAR  CYLINDER  MOOELS 

Circular  cylinders  with  diameters  ranging  in  size  from  2.5  in.  to 
5  in.  were  used  in  the  experiments.  The  cylinders  were  turned  on  a  lathe 


64 


65 


"1 


66 


M 


Fig.  13b  Back  view  of  the  tunnel  with  the  side  windows 


from  aluminum  pipes  and  polished  to  a  mirror-shine  surface.  The  length 
of  each  cylinder  was  such  as  to  allow  a  1/32  in.  clearance  between  the 
tunnel  wall  and  each  end.  As  part  of  the  total  length,  smooth  circular 
flanges  6.0  in.  in  diameter  and  1.8  in.  thick  were  attached  to  the  ends 
of  the  cylinder  to  eliminate  the  end  effects. 

Cylinders  were  hollow,  and  each  contained  a  beam  type  accelerometer. 
Figure  15  shows  a  cutaway  drawing  of  a  cylinder  with  the  accelerometer 
installed  in  the  exact  center  of  the  cylinder.  The  electrical  connection 
to  the  accelerometer  was  led. through  a  notch  in  the  male  fitting  at  the 
end  of  the  cylinder.  Sealing  of  the  unit  was  accomplished  by  0-rings 
and  threaded  fittings. 

The  same  cylinders  were  also  used  as  rough  cylinders  by  coating  them 
with  sand  of  desired  size.  The  relative  roughness  used  in  the  present 
experiments  was  kc  =  1/100.  For  this  purpose  clean  sand  was  screened 
through  the  use  of  a  combination  of  appropriate  size  of  sieves  to  separate 
the  desired  size.  The  cylinder  surface  was  then  cleaned  with  alcohol  and 
carefully  coated  with  an  epoxy  paint.  It  was  possible  to  obtain  an 
extremely  thin  and  uniform  coating  which  was  free  of  ridges,  waves  and 
local  buildup.  Next,  the  selected  sand  was  sprinkled  over  the  cylinder. 
The  excess  sand  was  easily  removed  after  15.  minutes  through  gentle  rubbing 
by  hand.  The  resulting  roughness  was  uniform  as  evidenced  by  the  two 
sample  photographs  shown  in  Figure  16.  The  pictures  of  the  rough  suface 
were  taken  with  a  scanning  electron  microscope.  The  actual  size  of  the 
sample  sand  is  shown  in  the  figure. 

C.  FORCE  MEASUREMENTS 

Two  identical  force  transducers,  one  at  each  end  of  the  cylinder, 
were  used  to  measure  the  instantaneous  force.  The  basic  transducer  was 
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Fig.  15 


manufactured  by  B.L.H.  Electronics,  Inc.,  under  the  trade  name  of  LBP-1. 
Two  sets  of  transducers  were  used  in  the  experiments.  Each  gage  in  one 
of  the  sets  had  a  capacity  of  500  1b  and  in  the  other  set  a  capacity  of 
250  lb.  The  overload  capacity  of  the  gages  was  200  percent.  The  deflec¬ 
tion  of  the  gage  under  500  lb  load  was  0.01  inch.  A  photograph  of  the 
basic  transducer  is  shown  in  Fig.  17.  The  bellows  which  protected  the 
strain  gages  had  to  be  waterproofed  in  such  a  manner  that  they  would  not 
adversely  affect  the  operation  of  the  gage  when  subjected  to  about  20 
feet  water  pressure  at  temperatures  ranging  from  65°F  to  160°F.  For  this 
purpose  the  bellows  were  completely  filled  with  Dow  Corning  340-RTV 
coating  without  bringing  the  rubber  into  contact  with  air  during  the 
filling  operation.  After  filling,  the  ends  of  the  bellows  were  sealed 
air  tight  with  special  clamps.  The  transducers  then  were  mounted,  as 
shown  in  Fig.  7,  to  the  aluminum  split-bearing  housing  blocks. 

It  is  seen  from  Fig.  18  that  the  transverse  force  acting  on  the 
cylinder  is  given  by 


Fl  =  F1  +  F2  -  MX  (2) 

in  which  MX  represents  the  inertial  force  due  to  the  acceleration  of 
the  cylinder.  The  sensitivity  of  the  accelerometer  was  adjusted  to  yield 
FL  =  0  while  the  cylinder  freely  oscillated  in  air.  This  method  allowed 
the  measurement  of  the  net  fluid  force  with  an  error  less  than  1.5  percent. 

D.  CYLINDER  DISPLACEMENT  MEASUREMENTS 

To  measure  the  transverse  displacement  of  the  cylinder  during  its 
oscillations  a  linear-variable-displacement  transducer  (LVDT)  was  installed 
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between  each  bearing  block  and  housing.  Each  LVDT  was  energized  by  a 
six  volt  DC  regulated  power  supply.  The  LVDT's  provided  displacement 
from  the  equilibrium  position  and  also  helped  to  verify  the  simultaneous 
and  in-phase  oscillation  of  both  ends  of  the  cylinder  (see  Figs.  19a  through 
19c  for  sample  displacement,  force,  and  acceleration  traces. 

E.  ACCELERATION,  ELEVATION,  OR  VELOCITY  MEASUREMENTS 

It  is  because  of  the  extreme  importance  of  the  accurate  measurement 
of  the  instantaneous  value  of  these  quantities  that  they  are  discussed 
here  separately « 

It  should  be  noted  that  the  measurement  of  the  amplitude,  accelera¬ 
tion,  elevation,  or  the  velocity  is  a  matter  of  interpretation  of  the 
signal  received  from  the  appropriate  transducer  in  light  of  one  of  the 
following  expressions. 


U 


m 


2ttA 
T  » 


_  dU  = 
dt 


(^-i2 


(3) 


in  which  T  is  the  tunnel  period  and  is  a  constant.  In  these  experiments 
use  was  made  of  a  differential  pressure  transducer  with  pressure  taps 
located  symmetrically  on  the  two  legs  of  the  tunnel  at  an  elevation  50 
inches  below  the  mean  water  level,  i.e.,  H  =  38  in.  Applying  Bernoulli's 
equation  for  unsteady  flow  between  each  pressure  tap  and  the  instantaneous 
level  of  water,  it  is  easy  to  show  that  twice  the  amplitude  of  the  free 
surface  oscillation  (virtual  amplitude)  is  given  by 

\  ■  2Ai  *  ■ AP  —  <4> 
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in  which  g  and  T  are  constant  and  H  is  kept  constant  (see  Fig.  14). 

Thus  the  signal  of  this  transducer  yielded  the  virtual  amplitude  or 
maximum  velocity  in  each  cycle.  It  was  entirely  free  from  noise  or 
free  surface  effects.  The  transducer  was  calibrated  and  its  linearity 
checked  before  each  series  of  experiments. 

F.  MEASUREMENT  OF  OTHER  BASIC  PARAMETERS 

The  total  stiffness  of  the  system,  Ks,  was  determined  by  suspending 
a  known  weight  from  the  center  of  the  cylinder,  and  the  displacement 
was  measured.  Then  K$  was  calculated  from  the  relation 


K  -  Force 
*s  Displacement 


(5) 


The  natural  frequency  in  air,  f  ,  and  the  natural  frequency  in 

ria 

water,  f  ,  were  determined  by  plucking  the  cylinder  and  allowing  it  to 
nw 

oscillate  freely  in  each  particular  fluid  medium.  The  displacement  was 
recorded  on  one  channel  of  a  three  channel  recorder.  A  sample  trace  is 
shown  in  Fig.  19d.  Along  with  the  frequencies,  the  logarithmic  decrement 
and  the  damping  factor,  r  ,  and  z,  ,  were  determined  in  air  and  water, 
respectively. 

Total  mass  of  the  system  was  measured  by  adding  up  the  mass  of  the 
components  of  the  oscillating  system.  The  system  mass  was  also  calculated 
from 
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to  check  the  correctness  of  the  frequency  evaluation.  The  results  were 
identical . 

Added  mass  M,  of  the  system  was  determined  from 

a 


M  +  M 

a 


(2'fnw)Z 


(7) 


For  the  purpose  of  verification  of  the  measured  values,  M  .  for  the 

a 

cylinder  was  computed  through  the  use  of  the  observed  values  of  frequencies 
and  found  to  be  within  one  tenth  of  one  percent  of  that  given  by  the 
potential  flow  theory  and  independent  experiments. 


G.  INSTRUMENTATION  AND  ELECTRONIC  CIRCUITRY 

The  output  of  each  of  the  two  LVDT's,  after  amplification  by  a 
low-gain  amplifier,  was  displayed  on  an  eight-channel  recorder  as  in 
Fig.  20.  One  of  the  two  signals  was  additionally  displayed  on  a  three- 
channel  recorder. 

The  output  of  the  force  transducers  was  separately  amplified  by  a 
carrier  amplifier,  filtered  by  a  Krohn-Hite  electronic  filter  and  dis¬ 
played  on  an  eight-channel  recorder.  Then  the  two  signals  were  added 
together  by  an  electronic  summing  circuit,  amplified  again  by  a  low 
level  preamplifier  and  displayed  on  the  eight-channel  recorder. 

The  accelerometer  output  was  similarly  amplified  by  a  carrier 
amplifier  and  filtered  by  a  low  pass  filter.  This  signal  was  also 
displayed  on  the  eight-channel  recorder.  It  was  then  subtracted  from 
the  sum  of  the  force  signals  through  the  use  of  a  summing-differencing 
electronic  circuit.  The  resultant  signal,  representing  the  fluid  force 
on  the  cylinder,  F^,  was  again  amplified  by  a  low  level  preamplifier  and 
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displayed  on  both  three-  and  eight-channel  recorders.  Figure  21  is  a 
block  diagram  of  the  circuit  described  above. 

The  signal  from  the  differential  pressure  transducer  was  of  such 
quality  that  it  did  not  require  any  filtering.  As  shown  in  Fig.  22, 
it  was  amplified  by  a  carrier  amplifier  and  displayed  on  the  two  afore¬ 
mentioned  recorders. 

H.  PROCEDURES 

Five  cylinders  were  used  in  these  experiments.  Additionally  two 
springs,  manufactured  to  the  specification  desired,  were  used.  The 
diameter  of  the  springs  were  sufficiently  large  to  work  freely  over  the 
2.5  in.  diameter  rod.  The  coil  diameter  of  the  springs  was  0.255  inches 
and  were  made  of  carbon  steel  wire.  Table  3  is  a  condensed  list  of  the 
cylinders  used,  kind  of  force  measured,  surface  condition,  and  spring 
constant  employed  in  the  experiment.  The  number  of  (x)  in  each  box 
indicates  the  number  of  experiments  performed  under  that  particular 
conditions.  The  check  mark  under  the  spring-constant  box  indicates  the 
kind  of  spring  used. 


Table  3 

Diam. 

LIFT 

DRAG 

K$  (Lb/ft) 

in. 

Smooth 

Rough 

Smooth  Rough 

418 

392 

2. 

X 

/ 

2.5 

X 

X 

/ 

3. 

XXX 

X 

X  X 

/ 

4 

XX 

XX 

/ 

/ 

5 

XXX 

XX 

/ 
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Fig.  21  Force  and  accelerometer  electronic  circuit 
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Fig.  22  Differential  pressure  transducer  and  circuitry 


In  some  of  the  experiments  it  was  found  advantageous  to  increase 
the  weight  of  the  oscillating  system  by  mounting  rectangular  lead  blocks 
to  the  side  of  the  oscillating  bearing  blocks  away  from  the  cylinder 
mountings.  Beside  adding  mass  and  therefore  lowering  the  natural  fre¬ 
quency,  it  was  discovered  that  adding  the  weight  at  this  location  also 
served  to  distribute  the  load  better  on  the  bearings,  thereby  reducing 
friction  and  damping. 

Procedures  for  an  experiment  were  as  follows: 

1.  The  cylinder  was  aligned  with  the  tunnel  by  adjusting  the 
bearings.  The  sensitivities  of  the  force  transducers  were  matched  by 
hanging  a  known  weight  from  the  exact  center  of  the  cylinder.  At  the 

same  time  a  calibration  factor  was  determined  by  measuring  the  displacement 
and  comparing  it  with  the  output  of  the  LVDT. 

2.  The  accelerometer  was  calibrated  by  plucking  the  cylinder  and 
adjusting  the  amplifier  gain  until  the  fluid  force  in  air  was  nulled  to 
within  1.5  percent  deviation. 

3.  f  .  and  c  were  also  determined  by  plucking  the  cylinder  in  air 

Ha  a 

from  the  trace  of  cylinder  oscillation  on  the  three-channel  recorder 
(see  Fig.  19d). 

4.  The  weight  was  then  suspended  from  the  cylinder  center,  this 
time  for  the  measurement  of  a  calibration  factor  in  lb/mm  of  the  three- 
channel  -recorder  output  at  a  particular  gain  setting. 

5.  The  U- tunnel  and  cylinder-support  housing  were  again  filled  to 
the  top  of  the  test  section,  and  all  air  was  vented  from  the  housings. 

6.  The  cylinder  was  then  lowered  to  the  center  of  the  test  section 
by  simultaneously  lowering  both  sets  of  spring  yokes,  each  the  same 
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number  of  turns,  to  compensate  for  the  buoyant  force  on  the  displaced 
cyl inder. 

7.  The  cylinder  was  plucked,  and  f  and  c  were  calculated  from 

nw  w 

the  free-oscil lation  traces. 

8.  The  access  to  the  horizontal  test  section  above  the  cylinder, 
through  which  the  weights  were  lowered  to  the  cylinder,  was  closed. 

The  tunnel  was  filled  with  water  from  a  sump  tank  located  near  the 
apparatus  or  from  the  city  water. 

9.  The  differential -pressure  transducer  was  bled  of  any  entrained 
air  and  its  calibration  verified. 

10.  Variable  speed  motor  control,  driving  the  air  valve,  was  adjusted 
to  yield  a  period  exactly  equal  to  that  of  the  tunnel. 

11.  Butteryfly-valve  on  top  of  the  air  chamber  and  the  gate-valve 
on  top  of  the  air  control  box  were  fully  closed.  The  gate  valve  at  the 
outlet  of  the  fan  was  fully  opened,  and  the  butterfly-val ve  on  top  of 
the  other  leg  was  fully  opened  before  the  fan  was  turned  on.  With  this 
arrangement  the  largest  amplitude  of  the  tunnel  oscillation  was  obtained. 

12.  Several  tunnel  oscillations  were  allowed  to  pass  in  order  to 
achieve  a  steady  state  tunnel  oscillation. 

13.  Amplifier  gain  settings  or  calibration  factors  were  recorded 
on  the  chart  paper  as  they  were  changed  as  desired. 

14.  Several  cycles  of  oscillation  were  recorded  at  a  chart  speed 
of  2  and  10  mm/sec  to  ensure  that  a  steady  state  is  reached  for  that 
particular  amplitude  of  the  tunnel. 

15.  Several  runs  were  made  at  a  chart  speed  of  25  mm/sec. 

16.  The  position  of  the  gate  valves  on  the  air  control  box  was 
changed  to  achieve  a  new  amplitude  of  the  tunnel  oscillation. 
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17.  Steps  14  through  16  were  repeated. 

18.  At  the  end  of  each  experiment,  f  ,  c  ,  f  ,  and  z,  were  measured 

nw  w  mu  o 

to  verify  the  initial  values. 

In  measuring  the  forces  in  the  in-line  direction,  steps  1  through 
4  were  repeated.  The  force  transducers  were  turned  90  degrees  and  pinned 
to  measure  the  in-line  force.  Then  the  steps  5  through  17  were  followed 
as  before. 
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III.  METHOD  OF  DATA  ANALYSIS 


A.  INTRODUCTION 

The  recorded  data  consisted  of  the  instantaneous  values  of  the 
transverse  displacement  of  the  cylinder,  lift  force,  and  the  acceleration 
of  flow.  Figures  19a  through  19c  show  a  sample  of  the  traces  of  the  above 
variables  during  two  cycles  of  flow  oscillation. 

The  data  were  subjected  to  extensive  analysis  in  terms  of  the  peak 
amplitudes  of  the  force  and  cylinder  displacement,  their  root-mean-square 
values,  Fourier  analysis,  spectral  analysis,  Morison's  equation  [96], 
and  finally,  in  terms  of  Duhamel's  superposition  integral. 


B.  ANALYSIS  OF  THE  PEAK  AMPLITUDES 

The  parameters  D,  K  ,  f  ,  f  ,  c  ,  and  c  were  recorded  for  each 
5  no  nw  q  w 

run.  Then  the  peak  values  of  the  transverse  force  (lift  force),  F^M; 
peak  values  of  the  cylinder  displacement,  XM;  and  the  periods  Tf  and 
Tc  for  the  force  and  displacement,  respectively,  were  evaluated  over  each 
half  cycle  of  flow  oscillation.  Subsequently,  the  variables  A^,  U^, 

K,  X^D,  Um/fnD,  and  C^M  were  calculated  for  further  study.  The  coefficient 
CLM  is  the  peak  lift  coefficient  and  is  given  by 


L  D 


(8) 


where  p  is  the  density  of  the  fluid;  and  L,  the  length  of  the  cylinder. 
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C.  DIGITIZATION  OF  THE  DATA 


The  evaluation  of  the  root-mean-square  values  of  the  various 
parameters,  Fourier  analysis,  spectral  analysis,  Duhamel's  superposition 
integral,  and  the  use  of  Morison's  equation  required  the  careful  digiti¬ 
zation  of  the  data. 

In  order  to  maintain  a  proper  relationship  between  the  phase  and 
amplitude  of  a  given  variable  a  suitable  frequency  of  digitization  or 
digitization  interval  had  to  be  chosen.  Otherwise  aliasing  could  occur 
which  is  a  potential  source  of  error  in  an  analog  to  digital  data  conver¬ 
sion  [103],  The  maximum  frequency  resolvable  is  called  the  Nyquist 
frequency  and  is  defined  as 


fNY  =  2At  ^ 

where  At  is  the  interval  of  digitization.  The  frequencies  in  the  original 
data  above  this  cut  off  frequency  are  folded  back  into  the  frequency  range 
from  zero  to  f^y  in  an  accordion  pleated  fashion.  Therefore,  in  order 
to  choose  a  suitable  At  the  maximum  frequency  encountered  in  the  experi¬ 
ments  had  to  be  known.  The  observation  of  the  analog  records  showed 
that  the  highest  frequencies  encountered  were  in  the  range  of  4  to  6  Hz. 

In  general,  it  is  a  good  rule  to  select  f^y  to  be  one-and-a-half  or 
two  times  greater  than  the  maximum  anticipated  frequency.  Therefore,  a 
Nyquist  frequency  of  10  Hz  was  chosen  from  which  an  interval  of  digitiza¬ 
tion  of  0.05  seconds  resulted.  This  corresponded  to  a  1.25  mm  interval 
on  the  chart  at  its  speed  of  25  mm/sec.  Any  digitization  interval  below 
this  value  would  have  given  correct  results.  The  smallness  of  this 
interval  is  defined  by  the  economy  of  the  effort. 
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In  case  of  the  use  of  data  in  a  spectral  analysis  the  number  of 
digitized  data  had  to  be  N  which  is  defined  as 

N  =  2  x  2P  p  =  1,  2,  3,  4,  ... 

This  consideration  and  others  such  as  the  ease  of  tapering  led  to  an 
interval  of  digitization  of  1.073  mm,  equivalent  to  At  =  0.043  seconds. 
This  time  interval  is  lower  than  the  conservative  value  of  0.05  seconds. 

The  length  of  the  records  were  different  but  never  less  than  three 
cycles  of  tunnel  oscillation  which  is  about  18  seconds.  This  resulted 
in  a  resolution  of  1/18  =  0.056  Hz.  The  digitized  records  were  long 
enough  to  satisfy  all  the  assumptions  inherent  in  the  analysis. 

The  analog  traces  at  a  fixed  tunnel  amplitude  were  digitized  by 
means  of  a  Hewlett-Packard  digitizer  model  9874A  which  was  linked  to  a 
Hewlett-Packard  9845A  computer.  The  movement  of  a  cursor  by  hand  over 
the  trace  generated  digitized  coordinates  of  a  given  curve.  The  informa¬ 
tion  was  then  stored  either  in  files  on  tapes  for  use  with  the  programs 
written  for  the  HP  computer  or  punched  on  cards  to  be  used  with  an  I.B.M. 
360  computer. 

D.  EVALUATION  OF  THE  ROOT -MEAN -SQUARE  VALUES 

The  use  of  digitized  data  together  with  a  suitable  computer  program 

yielded  the  rms  values  of  the  lift  coefficient,  (C.  )  and  the  relative 

L  rms 

cylinder  displacement,  (X)rn)S/D. 

The  record  length  for  the  evaluation  of  the  rms  values  was  at  least 
three  cycles  of  tunnel  oscillation.  A  careful  examination  of  the  data 
has  shown  that  a  record  length  of  three  cycles  was  more  than  sufficient. 


E.  SPECTRAL  ANALYSIS 


As  noted  earlier,  the  nature  of  the  lift  force  and  cylinder-displace¬ 
ment  traces  were  quite  regular,  particularly  in  the  synchronization  range. 
The  quasi -harmonic  shape  of  the  said  traces  required  a  deterministic 
rather  than  random  analysis.  However,  a  knowledge  of  the  frequency  con¬ 
tent  of  the  oscillations  was  of  primary  importance  in  determining  the 
modulations  in  the  amplitude  and  frequency  of  the  oscillation  of  the 
lift  force  and  cylinder  displacement.  Evidently,  only  through  spectral 
analysis  can  one  decompose  the  force  and  displacement  into  their  funda¬ 
mental  and  subharmonics  and  draw  conclusions  regarding  the  energy  content 
of  each  harmonic. 

The  following  is  a  condensed  description  of  the  methodology. 

If  x(t)  is  a  periodic  function  of  time  t,  with  a  period  T  ,  then 
it  could  always  be  expressed  as  an  infinite  trigonometric  series  of  the 
form 


x(t)  -a  +  I  [a  cos(^)  +  bn  sin(^)] 
0  n=l  "  's  "  's 


(10) 


where 


=  /  S  x { t )  dt 


on 


S  0 


a 

n 

n>l 


$-  I  *(0  cos(^)  dt 

SO  S 


02) 
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(13) 


T 

j-  /  S  x ( t )  sin(^-)  dt 
S  0  s 

The  mathematical  conditions  for  the  convergence  of  Eq.  (10)  are  extremely 
general  and  cover  practically  every  conceivable  engineering  situation 
[104]. 

The  Fourier  series  is  a  useful  tool  for  determining  the  frequency 
content  of  a  time-varying  signal.  However,  the  Fourier  series  always 
require  a  periodic  time  function.  To  overcome  this  shortcoming,  Fourier 
evaluated  his  series  as  he  let  the  period  of  the  waveform  approach 
infinity.  Then  the  Fourier  series  given  by  Eq.  (10)  turns  into  a  Fourier 
integral  and  Fourier  coefficients  represented  by  Eqs.  (11),  (12)  and 
(13)  become  continuous  functions  of  frequency,  called  Fourier  transforms. 
Therefore,  a  Fourier  integral  may  be  regarded  as  the  formal  limit  of  a 
Fourier  series  as  the  period  tends  to  infinity.  It  is  customary  to 
write  the  Fourier  series  in  their  complex  form.  Then  one  introduces 
the  Fourier  transforms  defined  by 

oo 

F  (f)  =  /  x(t)  e"i27Tft  dt  (Forward  Fourier  transform)  (14) 

-00 

00 

x(t)  =  /  Fx(f)  el2lTft  df  (Inverse  Fourier  transform)  (15) 


where 

e±i  2-n-f t  _  CqS ( 2mf t)  ±  i  sin(2Trft)  , 
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is  known  as  the  Kernel  of  the  Fourier  transform.  Now  Fx(f)  contains 
the  amplitude  and  phase  information  at  every  frequency  present  in  x(t) 
without  demanding  that  x(t)  be  periodic.  If  the  time  history  x(t)  of 
a  sample  function  is  not  periodic  it  cannot  be  represented  by  a  discrete 
Fourier  series.  Furthermore,  for  a  stationary  process,  x(t)  goes  on 
forever  and  the  necessary  condition 


/  |x(t)|  dt  <  «  (16) 

.00 

is  not  satisfied.  Unless  special  precautions  are  taken,  it  is  not 
therefore  possible  to  evaluate  the  Fourier  transform  of  x(t)  in  order 
to  obtain  information  about  the  frequency  composition  of  the  process. 

This  difficulty  can  be  overcome  by  Fourier  analysing,  not  the  sample 
function  of  the  process  itself,  but  rather  its  autocorrelation  function 
Rx(t). 

The  logic  behind  this  approach  is  that  the  autocorrelation  function 
gives  information  about  the  frequencies  present  in  a  random  process 
indirectly.  If  the  zero  value  of  the  process  x(t)  is  normalized  (or 
adjusted)  so  that  the  mean  value  of  the  process  is  zero  and  it  has  no 
periodicity,  then  the  autocorrelation  function  satisfies  Eq.  (16)  and 
its  Fourier  transform  could  be  written  as 


Vf)  »  /  Rx(t)  e 


-i2Trft 


dx 


(17) 


The  Sx(f)  is  called  mean  square  spectral  density,  energy  density,  or 


power  spectrum.  The  last  name  originated  from  electrical  engineering 


± 
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because  it  is  a  measure  of  the  power  loss  in  a  resistor.  It  can  be 
shown  that  Sx(f)  is  a  real,  even  function  of  frequency,  f  and  it  is 
always  positive  [103,  105]. 

One  sided  spectral  density  is  defined  as 

Gx(f)  »  2  Sx(f)  (18) 


It  can  be  shown  that 


\  =  /  G  (f)  df 

o  x 


(19) 


2 

where  a is  the  variance  of  the  process  [106]. 
The  one  sided  cross  spectrum  is  found  from 


vf)  ■  y  vt)*'1*rft,h 


(20) 


or 


Gxy<f)  =  Cxy<f>  •  1  Vf)  (21) 

where  C  (f)  ’s  the  co-spectrum.  It  is  the  real  and  in  phase  component. 

*y 

The  Q  (f)  is  the  quad-spectrum  and  forms  the  imaginary  part  or  out  of 

phase  component  of  G  (f). 

xy 

From  Eq.  (21),  the  phase  difference  between  the  process  x  and  y  is 


found  to  be 


(22) 


a. 


xy 


Q  ( f ) 

xyv  ' 

cm 

xy v  ' 


If  a  finite  record  of  length  T$  is  used,  it  can  be  shown  [103]  that  an 
estimate  of  the  cross-spectral  density  is  given  by 


(Ms) 


2  V(n)  X*(n) 

^  Af 


(23) 


where  X  and  Y  represent  the  finite  Fourier  transforms  of  x  and  y  respectively 
and  X*  is  the  complex  conjugate  of  X.  Then  an  estimate  for  spectral 
density  of  a  stationary  process  of  finite  length  would  be 


Gx(f,Ts) 


(24) 


where 


X(n)  =  l(an  -  ibn)  (25) 

in  which  an  and  bn  are  the  Fourier  coefficients  given  by  Eqs.  (11), 

(12),  and  (13).  If  N  digitized  data  are  used  to  evalute  an  and  bn 

2 

numerically  then  the  computation  time  is  proportional  to  N  .  Even  with 
high  speed  computers,  computation  of  discrete  Fourier  transforms  requires 
excessive  machine  time  for  large  N.  In  1965  Cooley  and  Tukey  published 
their  mathematical  algorithm  [107]  which  has  become  known  as  the  "fast 
Fourier  transform".  The  fast  Fourier  transform  (F.F.T.)  is  a  computational 
algorithm  which  reduces  the  computing  time  of  ap  and  bn  to  a  time 
proportional  to  N  log  N.  The  idea  of  F.F.T.  was  first  discovered  by 
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Runge  [108,  109]  in  1903  who  described  essentially  the  computation  of 

the  F.F.T.  algorithm  as  we  know  it  today. 

Any  spectral  estimate  has  only  two  degrees  of  freedom  and  the 

confidence  interval  turns  out  to  be  so  large  that  it  does  not  have  any 

statistical  meaning.  Therefore,  one  deduces  that  the  evaluation  of 

only  one  G  (f,T  )  for  a  record  length  does  not  have  any  significance. 

X  s 

There  are  several  methods  to  increase  the  number  of  degrees  of  freedom. 

A  simple  method  is  to  create  an  ensemble  of  records  each  of  which  has  a 
record  length  T  .  If  the  length  of  the  original  record  is  long  enough 
one  can  divide  the  record  into  k  records  each  T$  seconds  long.  Now 
ensemble  averaging  improves  the  degrees  of  freedom,  which  in  turn  improves 
confidence,  or  spectral  reliability  of  the  spectral  estimate  then 

G(f,T  )  =  1  l  G(f,T  k)  (26) 

S  K  k=1  s 

where  k  is  the  number  of  sample  function  forming  the  ensemble  and  the 
degrees  of  freedom  for  each  spectral  estimate  is  now  X  =  2k.  The  true 
spectrum  is  found  as 


G(f)  =  Lim  ~  l  G(f,T  ,k)  (27) 

T  -HO  K  k=l  s 

k-K» 

The  use  of  this  method  requires  that  the  process  be  ergodic.  This  means 
that  the  process  must  be  stationary,  i.e.,  for  any  sample  function  k, 
time  average  of  the  mean,  autocorrelation  and  all  other  statistical 
averages  do  not  differ  when  computed  over  different  sample  functions. 
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As  it  was  mentioned  earlier,  the  observation  of  the  force  and  cylinder 


displacement  traces  revealed  that  the  phenomenon  repeated  itself  at 
almost  every  three  cycles  of  flow  oscillation  for  a  given  Keulegan- 
Carpenter  number.  This  justified  the  assumption  of  ergodicity  for  the 
process. 

Two  or  more  sample  functions  each  of  which  comprising  at  least 
three  cycles  of  tunnel  oscillations  were  used  to  evaluate  spectral 
density  of  the  process  at  a  constant  Keulegan-Carpenter  number. 

A  flow  chart  of  the  program  used  to  evaluate  the  spectral  density 
is  given  in  Fig.  23.  The  data  was  first  "detrended"  to  remove  any  linear 
trenc  from  the  process.  Then  a  cosine  taper  window  was  used  to  reduce  the 
leakage  and  smoothen  discontinuities  at  the  "joints"  at  the  two  ends  of 
the  records. 

In  addition  to  the  spectral  density  for  each  data  set,  other  statisti¬ 
cal  averages  such  as  variance,  standard  deviation,  skewness,  and  kurtosis 
were  estimated.  Integration  of  spectral  density,  using  Eq.  (19),  was 
performed  to  compare  the  result  against  variance  which  was  evaluated 
directly  from  digitized  data.  The  closeness  of  these  two  was  also  a 
good  check  for  the  goodness  of  spectral  estimate. 

F.  IN  LINE  FORCE 

The  drag  coefficient  C.  and  the  inertia  coefficient  C  for  the  in- 
line  force  were  evaluated  in  order  to  determine  the  effect  of  the  cross- 
flow  oscillation  of  the  cylinder  on  the  in-line  force.  Data  reduction 
for  the  forces  in-line  with  the  direction  of  the  ambient  flow  was  based 
on  Morison's  equation  [96]  and  the  Fourier  analysis.  The  in-line  force, 
consisting  of  the  drag  force  F^  and  the  inertia  force  F^ ,  is  assumed 
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Fig.  23  Flow  chart  of  the  computer  program  for  the 
evaluation  of  the  spectral  density. 


n 


96 


to  be  given  by 


F  ’  Fd*F<  '  ?CdL°p  |U|  (28) 

in  which  Cd  and  Cm  represent  respectively  the  drag  and  inertia  coeffi¬ 
cients  and  U  is  the  instantaneous  velocity  of  the  ambient  flow.  For  an 
harmonic  flow  represented  by 

U  = 

where 

e  = 

equation  (28)  may  be  written  as 

-  Cd  | cos  8 (  cos  0  +  Cm  y  sin  0  (30) 

and  Cm  were  evaluated  from  (30)  by  Keulegan 

,  2ir  F 

*  -  j  /  - 5 - cos  0  d0  (31) 

■  T~  /2” - - sin  6  de  (32) 

0  »V01 

measured  force. 


T - 2 - 

b  VDL 

The  Fourier  averages  of 
and  Carpenter  as  (111] 


in  which  F  represents  the 


-  IL  cos  0  (29) 
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The  In-line  and  lift  forces  were  measured  by  two  force  transducers 
at  each  end  of  the  cylinder  while  It  was  oscillating  in  the  transverse 
direction.  The  three  channel  recorder  with  three  strip  chart  paper  was 
used  to  record  forces  and  the  flow  oscillation.  Figures  24  and  25  show 
sample  in-line  and  transverse  force  traces.  The  in-line  force  was  then 
digitized  every  0.1  seconds  and  the  digitized  data,  along  with  other 
parameters  such  as  0,  L,  v  ,  and  Am,  were  fed  into  a  computer  program  to 
evaluate  and  Cm  from  Eqs.  (31)  and  (32). 

G.  DUHAMEL'S  SUPERPOSITION  INTEGRAL 


One  of  the  methods  used  to  investigate  the  relationship  between  the 
instantaneous  force  and  displacement  in  a  linear,  oscillating,  system 
is  the  so-called  Duhamel's  superposition  integral.  Assume  that  the 
oscillation  of  the  cylinder  is  represented  by 


MX+2MswnX+ksX  =  FL(t)  (33) 


where  F^(t)  is  the  instantaneous  lift  force.  One  can  imagine  that  the 
force  FL ( t )  is  composed  of  a  series  of  impulses  [111].  Then  through 
the  use  of  the  convolution  or  superposition  integral  one  can  evaluate 
the  instantaneous  response  of  the  system  to  the  force  FL(t)  from  [111] 

t  -?<*>n(t-tv) 

-  (Xft  -  j jL-  /  - sinfl-?  wntvdtv)  cos  o,nt 


n  o 


Vi  -  c2 


*  *X  ‘  t  -?wn(t-t  ) 

♦  (iu  i  fat 


’w  Mui 
n  n  o 


cos  y^7  untvdtv)s1nu>nt  (34) 
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acceleration 


Fig.  25  A  sample  of  the  transverse  force  trace 


In  which  XQ  and  XQ  denote  the  Initial  condition  and  ty  Is  the  variable 
time. 

For  this  purpose,  the  lift  force  trace  was  digitized  every  .02 
seconds  and  used  as  the  Input  to  Eq.  (34)  along  with  other  parameters 
characterizing  the  oscillating  system  such  as  c,  0,  M,  etc.  A  computer 
program  then  performed  the  calculations.  This  computation  served  also 
as  a  confirmation  for  the  consistency  of  the  lift  force  data. 


IV.  RESULTS 


A.  INTRODUCTION 

The  experimental  results  will  now  be  reported,  showing  first  the 
dependence  of  the  relative  displacement  and  lift  coefficient  on  the 
flow  parameters,  and,  secondly,  the  analysis  of  the  data  through  the  use 
of  the  Fourier  analysis,  spectral  analysis,  and  Duhamel's  integral. 

Lastly,  the  in-line  force  acting  on  transversely  oscillating  cylinders 
will  be  presented. 

B.  REPRESENTATIVE  RELATIVE  DISPLACEMENTS  AND  LIFT  COEFFICIENTS 

The  variation  of  the  relative  displacement  XM/D  as  a  function  of 
K  -  Um  T/D,  Ur  =  Um/fnwD,  and  fc/fnw  is  shown  in  Figs.  26  through 

29  for  a  3  in.  smooth  cylinder.  Data  for  other  representative  smooth 
and  rough  cylinders  are  shown  in  Figs.  30  through  45  in  Appendix  A.  In 
the  foregoing,  Um  represents  the  maximum  velocity  of  flow  during  a  cycle; 
f^,  the  frequency  of  the  transverse-force  oscillations  between  two 
successive  peaks;  f  ,  the  natural  frequency  of  the  cylinder  in  still 
water;  and  f  ,  the  frequency  of  the  cylinder  oscillations.  The  relative 
roughness  for  all  rough  cylinders  was  kc/D  *  1/100  where  kc  is  the 
average  sand- roughness  height. 

The  variation  of  the  rms  values  of  X/D  with  K,  Ur,  and  ^</fnw 

for  the  3  in.  smooth  cylinder  is  shown  in  Figs.  46  through  49.  Similar 
data  for  other  cylinders  are  shown  in  Figs.  50  through  65  in  Appendix  B. 

The  peak  amplitudes  of  the  lift  coefficient  are  shown  in  Figs.  66  through 
69  as  a  function  of  K,  U  ,  fj f-..  and  ff/f  for  the  3  in.  smooth  cylinder. 

r  T  nw  i  w 
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Data  for  other  cylinders  are  shown  In  Figs.  70  through  85  In  Appendix  C. 

The  rms  values  of  for  the  3  In.  cylinder  are  shown  In  Figs.  86  through 
89  for  the  3  In.  cylinder.  The  remainder  Is  shown  In  Figs.  90  through 
105  In  Appendix  0. 

The  data  presented  in  the  figures  cited  above  exhibit  certain 
characteristics  which  will  be  discussed  below.  Firstly,  all  the  plots 
showing  the  peak-amplitude  values  of  the  relative  amplitude  and  the  lift 
coefficient  show  considerable  scatter.  This  is  partly  due  to  the  amplitude 
and  phase  modulations  in  the  force  and  displacement  traces  and  partly  due 
to  the  inherently  random  nature  of  the  vortex  shedding  in  steady  and 
harmonic  flow  about  stationary  and  oscillating  cylinders.  The  rms  values 
of  the  displacement  and  force  exhibit  much  smoother  variation  with  respect 
to  the  parameters  K,  U  ,  f-/f_UJ«  and  f./f  .  The  data  span  a  broad  range 
of  Keulegan-Carpenter  numbers.  Within  this  range  the  relative  amplitude 
shows  one  dominant  and,  at  times,  one  secondary  peak.  For  example.  Fig. 

26  shows  that  XM/D  reaches  a  maximum  value  of  0.27  at  a  Keulegan-Carpenter 
number  of  about  55.  A  second  and  considerably  smaller  peak  occurs  at 
about  K  =  135  (see  Figs.  26  and  46).  The  region  in  which  the  first  peak 
occurs  is  called  the  primary  region  of  synchronization  and  the  region  in 
which  the  secondary  peak  occurs  is  called  secondary  synchronization. 

A  tertiary  synchronization  region  observed  in  steady  flow  past  oscillating 
cylinders  has  not  been  found  in  the  harmonically  oscillating  flow.  Clearly 
the  primary  synchronization  region  extends  over  a  small  range  of  K  values. 
Outside  the  synchronization  region  the  cylinder  continues  to  oscillate 
at  small  amplitudes  as  a  body  subjected  to  forced  oscillations  by  shedding 
of  the  vortices.  The  dependence  of  the  primary  synchronization  region  on 
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the  reduced  velocity  Um/fnwD  is  seen  clearly  in  figures  similar  to 
Figs.  27  and  47.  The  primary  synchronization  (referred  to  hereafter  as 
the  synchronization)  occurs  at  Ur  values  characteristic  of  those 
encountered  in  steady  flows  (see  Fig.  3  where  V/ ( 2-rrf  0)  =  1].  The 

nw 

frequency  ratios  of  fJ f  and  f  /f  (see  Figs.  28  and  29  and  Figs.  48 

T  I  ItI  ( IW 

and  49)  show  that  synchronization  takes  place  when  the  frequency  of  force 
oscillations  nearly  coincides  with  the  natural  frequency  of  the  cylinder. 
Furthermore,  the  cylinder  oscillates  at  a  frequency  close  to  its  natural 
frequency. 

The  lift  coefficient  at  synchronization  is  several  times  larger 
than  that  in  the  non-synchronization  region  (see,  for  example.  Figs. 

66-69  and  Figs.  86-89).  As  will  be  noted  later,  the  said  lift  coefficient 
is  2.75  times  that  for  the  stationary  cylinder  at  the  corresponding 
Keulegan-Carpenter  number.  The  amplification  of  the  lift  force  is  a 
consequence  of  the  non-linear  interaction  between  the  flow  and  the  cylinder 
as  discussed  in  Chapter  I. 

As  noted  earlier,  the  random  nature  of  the  vortex  shedding  gives 
rise  to  modulations  in  both  displacement  and  force  traces.  It  is, 
therefore,  useful  to  examine  the  relative  influence  of  the  harmonics  of 
the  displacement  and  lift  force. 

C.  FOURIER  ANALYSIS  OF  THE  SYNCHRONIZATION  REGION 

When  a  signal  is  periodic  over  T$  it  is  convenient  to  use  a  Fourier 
series  representation.  The  dependence  of  the  Fourier  coefficients  on 
the  pertinent  dimensionless  groups  may  then  be  used  to  describe  the 
signal  completely.  Even  when  the  signal  is  not  completely  regular,  it 
is  still  useful  to  investigate  the  contributions  to  it,  associated  with 


various  harmonics  of  the  wave  frequency.  Thus,  if  the  displacement  and 
force  signals  are  considered  as  given  approximately  by 

oo 

X(t)  =  l  Xn  cos(27rnft/T  +  <f>  )  (35) 

n=l  n  n 

and 

00 

C.  (t)  =  l  C.  cos(2Trnft/T)  (36) 

L  n=l  Ln 

in  which  Xn  and  C^n  are  the  effective  Fourier  components  and  may  be 
deduced  through  the  use  of  a  suitable  computer  subroutine. 

Figures  106  through  113  show  the  normalized  displacement  and  their 
corresponding  normalized  Fourier  coefficients  for  all  smooth  and  rough 
cylinders.  An  examination  of  the  Fourier-coefficient  plot  shows  that  the 
predominant  frequency  invariably  occurs  at  f  =  9f  -  9/T  where  T  is  the 

W 

period  of  flow  oscillation  in  the  tunnel.  The  experiments  by  Sawaragi 
[99]  have  also  shown  that  the  predominant  frequency  is  equal  to  a  multiple 
of  the  flow  oscillation  frequency.  The  relationship  between  the  predominant 
frequency  and  Strouhal  number  will  be  discussed  shortly.  Figures  106 
through  113  also  show  that  there  are  secondary  harmonics  which  are 
considerably  smaller  than  the  predominant  harmonic.  A  correct  represen¬ 
tation  of  the  phase  and  amplitude  modulations  of  the  displacement  requires 
the  use  o*  at  least  two  or  three  harmonics.  Figures  114  through  121  show 
representative  normalized  force  traces  and  the  corresponding  normalized 
Fourier  coefficients. 

Once  again  it  is  found  that  the  predominant  frequency  of  the  lift 
force  is  ff  =  9  fw.  Also,  there  are  relatively  strong  higher  harmonics 


121 


Fig.  106  Normalized  displacement  and  Fourier  coefficients 
for  3  in.  smooth  cylinder,  c  =  0.052. 
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Fig.  107  Normalized  displacement  and  Fourier  coefficients 
for  3  in.  rough  cylinder,  e  =  0.067. 


Fig.  108  Normalized  displacement  and  Fourier  coefficients 
for  4  in.  smooth  cylinder,  c  ■  0.046. 
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Fig.  109  Normalized  displacement  and  Fourier  coefficients 
for  4  In.  rough  cylinder,  ?  =  0.030. 
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Fig.  112  Normalized  displacement  and  Fourier  coefficients 
for  5  In.  smooth  cylinder,  C  *  0.055. 
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Fig.  113  Normalized  displacement  and  Fourier  coefficients 

for  5  in.  smooth  cylinder,  c  *  0.054,  (with  lead  blocks) 
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Fig.  114  Normalized  lift  force  and  Fourier  coefficients 
for  3  in.  smooth  cylinder,  5  *  0.052. 
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Fig.  116  Normalized  lift  force  and  Fourier  coefficients 
for  4  In.  smooth  cylinder,  5  =  0.046. 
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Fig.  117  Normalized  lift  force  and  Fourier  coefficients 
for  4  in.  rough  cylinder,  x,  =  0.060. 
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Fig.  120  Normalized  lift  force  and  Fourier  coefficients  for 

5  in.  smooth  cylinder,  5  =  0.055,  (with  lead  blocks). 
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Fig.  121  Normalized  lift  force  and  Fourier  coefficients  for 
5  in.  smooth  cylinder,  c  =  0.054,  (with  lead  blocks) 
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on  either  side  of  the  predominant  frequency.  It  is  found  that  the 

instantaneous  lift  force,  as  well  as  displacement,  may  be  represented 

fairly  well  through  the  use  of  their  predominant  and  two  or  three 

additional,  most  significant,  harmonics.  Figures  122a  through  122c  show 

the  Fourier  representation  of  a  sample  displacement  trace  for  a  4  in. 

rough  cylinder  with  one,  two,  and  three  Fourier  components,  respectively. 

Similar  results  for  the  force  trace  are  shown  in  Figs.  123a  -  123e. 

Evidently,  the  displacement  is  fairly  well  represented  by  the  first 

three  components.  The  force  trace  for  the  particular  cylinder  required 

the  use  of  four  Fourier  components. 

The  occurrence  of  the  predominant  frequency  at  9f  is  not  without 

w 

a  reason.  If  the  Strouhal  number  is  defined  in  terms  of  the  maximum 
velocity  Um,  then  one  has 

fv° 

S  *  (37) 

m 

which  varies  from  about  0.1  to  0.23  during  a  cycle  in  harmonic  flow. 
Multiplying  the  numerator  and  the  denominator  by  T,  the  period  of  the 
flow  oscillation,  one  has 


fvDT  .  fvT 
‘  K 


(38) 


It  is  fully  realized  that  the  vortex  shedding  frequency  f  does  not 
remain  constant  throughout  the  cycle.  In  fact,  f  is  largest  near  the 
instant  of  Um  and  decreases  as  U  approaches  zero.  Consequently, the 
fundamental  question  as  to  which  representative  vortex  shedding  frequency 
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the  cylinder  oscillation  synchronizes  may  be  answered  only  through  the 
Fourier  analysis.  Replacing  f  by  N/T  one  has  SK  =  N.  In  the  present 
investigation  the  synchronization  occurs  in  the  range  of  Kelulegan-Carpenter 
numbers  from  about  48  to  60.  With  N  =  9,  one  has  S  =  0.15  -  0.18. 

Using  a  value  of  S  =  0.16  for  the  perfect  synchronization  (at  which  the 
displacement  is  maximum)  one  finds  K  =  56.  In  other  words,  the  cylinder 
undergoes  synchronized  oscillations  when  the  natural  frequency  of  the 
cylinder  coincides  with  the  vortex  shedding  frequency  based  on  a  Strouhal 
number  of  0.16,  as  verified  by  the  experiments.  Consequently  the  integer 
N  may  be  smaller  or  larger  depending  on  the  Keulegan-Carpenter  number 
at  which  synchronization  is  realized.  For  example,  for  a  cylinder  whose 
natural  frequency  is  twice  the  flow  oscillation  frequency  (N  =  2),  the 
synchronization  takes  place  at  K  =  12.5.  This  result  is  in  conformity 
with  the  results  of  Sawaragi  [99].  However,  the  Keulegan-Carpenter  number 
in  his  experiments  is  based  on  the  rms  value  of  the  wave  induced  veloci¬ 
ties  along  the  vertical  cylinder.  Thus,  an  exact  verification  of  the 
corresponding  results  is  not  possible.  Finally,  it  should  be  noted 
that  the  use  of  the  rms  value  of  the  flow  velocity  and  the  frequency  at 
which  perfect  synchronization  takes  place  results  in  S  =  0.22. 

D.  SPECTRAL  ANALYSIS 

The  spectral  density,  coherence,  and  phase  at  perfect  synchronization 
are  shown  in  Figs.  124  and  125  for  the  3  in.  smooth  and  4  in.  rough 
cylinders.  Similar  results  in  the  non-synchronization  region  are  shown 
in  Figs.  126  and  127  for  a  representative  value  of  K  =  80.  The  spectral 
density  for  a  broad  range  of  K  values  for  the  3  in.  cylinder  is  shown  in 
Figs.  128  through  143  in  Appendix  E.  Similar  data  for  the  4  in.  rough 
cylinder  are  shown  in  Figs.  144  through  158  in  Appendix  F. 
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PHASE  COHERENCE  LOG  SPECTRAL  DENSITY 


Fig.  124  Spectral  analysis  for  a  3  in.  smooth  cylinder 
at  resonance,  e  =  0.052,  K  =  51. 
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,  It.is  c1ea«i»r»'om  F.igs.  ^ 24-  and  125-that  the  spectra’  aenstty'is  *® 
sharply  defined  at  synchronization  and  the  force  and  displacement  peaks 
occur  at  the  same  frequency.  Furthermore,  the  coherence  is  nearly  per¬ 
fect.  In  the  non-synchronized  region  the  spectral  density  for  the  force 
becomes  broad  band  and  the  peaks  of  the  displacement  and  force  do  not 
occur  at  the  same  frequency  as  seen  from  Figs.  126  and  127.  The  peak 
of  the  displacement  is  relatively  better  defined  and  continues  to 
occur  near  the  natural  frequency  of  the  cylinder  (f  =  1.69  for  the  3 
in.  smooth  cylinder  and  f  =  1.54  for  the  4  in.  rough  cylinder).  The 

nw 

results  of  the  spectral  analysis  confirm  the  findings  of  the  Fourier 
analysis  that  at  resonance  the  force  and  displacement  maximums  occur  at 
the  same  frequency. 

E.  DUHAMEL’S  INTEGRAL 

Figures  159  through  163  show  a  comparison  of  the  measured  and  calcu¬ 
lated  values  of  the  instantaneous  displacement  as  a  function  of  time. 

The  measured  force  shown  in  the  upper  half  of  each  figure  was  used  as 
input  in  Duhamel 's  integral,  as  given  by  Eq.  (34).  The  close  agreement 
between  the  measured  and  calculated  displacements  is  an  indication  of  the 
proper  measurement  of  the  various  quantities  such  as  mass,  damping  ratio, 
frequency,  etc. 

A  sensitivity  analysis  was  carried  out  by  varying  the  frequency  f 
in  Eq.  (34)  while  maintaining  all  other  input  parameters  unaltered.  The 
purpose  of  this  analysis  was  partly  to  explore  whether  the  natural  fre¬ 
quency  of  the  oscillating  cylinder  remained  equal  to  that  determined  in 
still  water  and  whether  the  prediction  of  Duhamel's  integral  strongly 
depended  on  the  natural  frequency.  The  implication  of  the  equivalence 


of*the-rofctfral  ffCquCTcy  of  tfTfe  cyTindeT  undergoing  transverse  oscilla¬ 
tions  in  harmonic  flow  to  that  in  still  fluid  is  that  the  added  mass  of 


the  synchronized  cylinder  is  identical  to  the  displaced  mass  of  the 
cylinder.  If,  on  the  other  hand,  the  natural  frequency  of  the  synchronized 
cylinder  were  to  be  higher  than  that  in  still  water,  one  would  conclude 
that  the  added  mass  of  the  synchronized  cylinder  is  smaller  than  its 
displaced  mass.  Extensive  calculations  with  various  values  of  the 
natural  frequency  have  shown  that  the  natural  frequency  of  the  synchronized 
cylinder  is  5  percent  larger  than  that  for  the  natural  frequency  in 
still  water.  The  difference  between  the  two  frequencies  may  be  partly 
due  to  the  experimental  errors  and  partly  due  to  the  difference  between 
the  actual  added  mass  and  displaced  mass.  King  [112]  conducted  a  series 
of  experiments  in  still  water  (by  plucking  a  cantilevered  cylinder)  and 
also  in  flowing  water  at  the  fundamental  mode  of  the  flow-excited  vibra¬ 
tions.  From  a  comparison  of  the  measured  frequencies  he  concluded  that 
the  still-water  and  flow-excited  frequencies  are  virtually  equal,  i.e., 
the  added  mass  is  unaffected  by  streaming  flow.  Apparently,  this  conclu¬ 
sion  is  valid  only  at  the  resonance  condition.  The  results  of  the  present 
investigation  also  showed  that  the  added  mass  is  relatively  unaffected 
by  the  harmonically  oscillating  flow.  One  must,  however,  note  that  the 
best  agreement  between  the  measured  and  calculated  displacement  is 
obtained  with  a  natural  frequency  approximately  5  percent  larger  than 
the  still -water  natural  frequency. 

The  foregoing  also  points  out  the  important  fact  that  a  simple  com¬ 
parison  of  the  natural  frequencies  measured  in  still  water  and  in  flowing 
water  at  resonance  is  not  sufficient  for  an  accurate  assessment  of  the 
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Fig.  159  Comparison  o*  measured  and  calculated  displacement  by 

Ouhamel's  Integral  for  a  3  in.  smooth  cylinder,  ;  *  0.052. 
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and  calculated  displacement  by 
a  4  in.  smooth  cylinder,  c  =  0.046. 


ague 


HB 


Fig.  161  Comparison  of  measured  and  calculated  displacement  by 

Duhamel's  Integral  for  a  4in.  rough  cylinder,  c  =  0.030. 
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Fig.  162  Comparsion  of  measured  and  calculated  displacement  by 

Duhamel's  integral  for  a  4in.  rough  cylinder,  5  =  0.060. 


and  calculated  displacement  by 
a  5in.  smooth  cylinder,  c  *  0.054. 


-  vacation  or  ]«!;  of^ariation  of  theadded  mass.  The  use* of  Duhamers 
integral  is  a  more  reliable  method  for  the  determination  of  the  natural 
frequency  at  resonance.  Be  that  as  it  may,  it  is  rather  surprising  that 
the  added  mass  of  a  body  oscillating,  at  its  natural  frequency,  in  a 
fluid  otherwise  at  rest,  is  nearly  the  same  with  that  undergoing  syn¬ 
chronized  oscillations  in  a  fluid  in  motion. 

F.  DRAG  COEFFICIENT  FOR  TRANSVERSELY  OSCILLATING  CYLINDERS 

Experiments  with  cylinders  undergoing  transverse  oscillations  in 
steady  flow  have  shown  that  [113]  the  in-line  force  increases  with  X/D 
since  the  cylinder  presents  a  larger  apparent  projected  area  to  the  mean 
flow. 

From  an  engineering  viewpoint,  the  significance  of  the  magnitude  of 
the  in-line  force  is  that  a  cylinder  or  cable  excited  by  the  flow  to 
oscillate  in  the  transverse  direction  may  be  subjected  to  in-line  forces 
larger  than  those  assumed  in  its  design.  Furthermore,  the  deflections 
caused  by  the  in-line  force  of  sufficiently  flexible  cylinders  tend  to 
couple  with  transverse  oscillations  and  affect  not  only  the  magnitude  of 
the  transverse  oscillations  but  also  the  path  of  the  cylinder  motion. 

It  is  because  of  these  reasons  that  the  in-line  force  acting  on  the 
cylinders  undergoing  oscillations  in  the  transverse  direction  to  the 
harmonic  flow  has  been  measured.  The  drag  coefficient  has  been  determined 
through  the  use  of  Eq.  (31)  and  is  shown  in  Figs.  164  through  166  for 
three  representative  cylinders.  Also  shown  in  these  figures  are  the  mean 
lines  for  the  drag  coefficient  of  the  corresponding  stationary  cylinders. 

In  general,  the  oscillating-cylinder  drag  coefficients  are  larger  than 
those  for  the  stationary  cylinder.  However  the  increase  is  not  as  dramatic 
as  in  the  case  of  steady  flow  about  vibrating  cylinders  [114]. 
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The  inertia  coefficient  has  j^lsQ.-beeji.pva.lijateji  an*LJt  atas  icur.d  — 
that  it  did  not  differ  from  that  for  stationary  cylinders.  This  is 
primarily  because  of  the  fact  that  at  Keulegan-Carpenter  numbers  larger 
than  about  30  the  in-line  force  is  primarily  drag  dominated  and  the 
effect  of  the  inertial  force  is  negligible  [see  Eq.  (30)]. 
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V.  ANALYSIS  AND  DISCUSSION  OF  RESULTS 


A.  INTRODUCTION 

In  the  following,  the  results  obtained  with  various  analyses  and 
models  will  be  compared  with  those  obtained  experimentally.  First  the 
general  equation  of  motion  will  be  presented  and  then  the  oscillator 
model  of  Hartlen  and  Currie  [87]  will  be  discussed  with  reference  to 
the  experimental  data  presented  in  the  previous  chapter. 

B.  EQUATION  OF  MOTION  AND  COMPARISON  OF  RESULTS 

The  equation  of  motion  for  an  elastically-mounted  and  linearly- 
damped  cylinder  may  be  written  as 

MX  +  2Mu)nsX  +  Ks  =  FL(t)  (39) 

where  X  represents  the  transverse  displacement  of  the  cylinder;  M,  the 
mass  of  the  cylinder  plus  the  displaced  mass;  the  damping  ratio  in 
still  water;  ton,  the  natural  frequency  of  the  synchronized  cylinder; 

K^,  the  spring  constant;  and  FL(t),  the  exciting  force.  FL(t)  may 
be  presented  by 

FL(t)  =  D  L  (C^  Um2  Clim\-)ne  (40) 

in  which  (C^#)  is  the  lift  coefficient  for  the  stationary  cylinder 
at  the  corresponding  Keulegan-Carpenter  number  and  Reynolds  number. 

The  ratio  C^UJ/fC^0)  represents  the  amplification  of  the  lift  force 


164 


r 


at  or  in  the  neighborhood  of  perfect  synchronization.  Even  during  the 
period  of  nearly  perfect  synchronization  the  exciting  force  does  not 
retain  a  constant  amplitude  and  frequency.  This  fact  has  already  been 
demonstrated  through  the  use  of  the  Fourier  analysis.  In  general,  it 
is  necessary  to  consider  the  harmonics  of  the  lift  force.  Rewriting 
Eq.  (36),  one  has 

00 

C.  (t)  =  l  C.  cos{2imfft/T)  (36  -  repeated) 

L  n=l  Ln  T 

Inserting  CL(t)  in  Eqs.  (39)  and  (40),  one  obtains  a  differential 
equation  which  may  be  solved  by  assuming 

oo 

X(t)  xn  cos(2irnfft/T+*  )  (35  -  repeated) 

n=l 

The  general  solution  of  Eq.  (39)  is  given  by 


jrI£^CL#)rms^UnD^2cLn^CL°^nns 

{[1  -  (— )2]2  +  (2C  naif)2)1/2 
The  relative  amplitude  X*n/D  is  thus  given  by 


cos(ntdft/T+<!>n) 


(41) 


r  ?'£TT<CL”)™S>(^r>2W<CL°>nnS 

ji  =  _ n _ 

D  rh)*  n  f>  —  '  nau  «  1/2 

([1  -  (-1)2]2  +  (2  c  r-1)  > 
wn  “n 


(42) 
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Finally,  the  phase  angle  reduces  to 


nto, 

2  ?  — 

tan  1  [—  n 


ntO,  p 

1  -  (— f-)Z 
'go  ' 
n 


(43) 


The  amplification  ratios  ^i_n/(^|_°)rms  ^ave  ^een  evaluated  for  the 
fundamental  and  two  harmonics  through  the  use  of  the  Fourier-decomposed 
exciting  force.  The  results  have  shown  that  the  said  ratio  for  the 
predominant  harmonics  at  synchronization  are  given  by 


CLp/<CL°>™s  *  2-75 

±  0.15 

(44a) 

■  '-15 

±  0.10 

(44b) 

V'VW  *  °-70 

±  0.07 

(44c) 

^CL*t**rms/^Cl”*ras  " 

2.80  ±  0.20 

(44d) 

in  which  pi  represents  the  predominant  harmonics  (see  Figs.  114-121). 
The  ( C|_ ° ) ^5  values  for  the  corresponding  stationary  cylinders  were 
taken  from  the  data  presented  in  1101]. 

It  is  rather  remarkable  that  the  amplification  ratios  at  synchroni¬ 
zation  remain  fairly  constant  in  spite  of  the  wide  range  of  the  test 
parameters.  Consequently,  they  can  be  used  with  confidence  in  the  pre¬ 
diction  of  the  relative  amplitude  of  the  cylinder  oscillations.  It 
should  be  noted,  in  passing,  that  there  is,  at  present,  no  theoretical 
means  to  predict  the  amplitude  ratios.  The  use  of  the  Navier-Stokes 


equations  is  limited  to  very  small  Reynolds  numbers.  The  discrete-vortex 
model  is  limited  to  steady  flow  about  stationary  and  oscillating  cylinders 
[10].  It  is  not  currently  possible  to  apply  the  discrete-vortex  model 
to  the  hydroelastic  oscillation  of  cylinders  in  harmonically  oscillating 
flow.  This  is  primarily  because  of  the  difficulty  of  the  determination 
of  the  numerous  and  mobile  separation  points  on  the  cylinder.  Consequently, 
one  must  rely  on  the  experimentally  determined  values  of  the  amplifi¬ 
cation  ratios  in  predicting  the  cylinder  response. 

Equations  (44a)  through  (44c)  together  with  Eq.  (41)  have  been  used 
to  predict  the  response  of  all  the  cylinders.  Figures  167  and  168  show 
two  of  these  plots  where  the  predicted  and  measured  responses  are  seen 
to  compare  reasonably  well.  Similar  comparisons  with  all  other  cylinders 
proved  to  be  equally  satisfactory. 

A  critical  examination  of  the  relative  significance  of  the  second 
and  third  harmonics  on  the  cylinder  response  has  shown  that  the  response 
is  largely  determined  by  the  predominant  harmonic.  This  has  already  been 
anticipated  because  of  the  fact  that  the  second  and  third  harmonics  of 
the  relative  displacement  are  an  order  of  magnitude  smaller  than  the 
predominant  harmonic  (see  Figs.  106-113).  Assuming  that  the  fundamental 
harmonic  of  the  exciting  force  is  representable  by  (C[_)nns  and  that 
*  “n  at  synchronization,  Eq.  (42)  reduces  to 


X 

D 


rms| 

16rT2Me;/[^LDZ(CLo)rms] 


(45) 


A  comparison  of  Eqs.  (44a)  and  (44d)  shows  that 


<CLW<CL°>nns  5  CLp/<CL°W  <46> 

Thus,  it  is  reasonable  to  expect  that  Eq.  (45)  will  yield  the  fundamental 
harmonic  of  the  relative  displacement  fairly  accurately.  As  noted 
earlier,  higher  order  harmonics  of  X/D  are  an  order  of  magnitude 
smaller  and  certainly  quite  negligible. 

Equation  (45)  shows  that  a  relationship  should  exist  between  X/D 
and  the  response  parameter 


R 


P 


.  .M-S _ 

L  <V>ms 


(47) 


for  both  smooth  and  rough  cylinders.  Figure  169  shows  the  relationship 
between  X/D  and  Rp  for  all  test  cylinders.  The  relationship  between  the 
two  parameters  is  rather  remarkable  in  view  of  the  fact  that  only  the 
fundamental  harmonic  of  the  exciting  force  has  been  used  and  it  was 
assumed  that  =  un  =  u  .  Figure  169  also  shows  that  the  response  of 
a  given  body  to  a  given  flow  must  necessarily  depend  on  the  starting 
dynamical  conditions,  i.e.,  on  the  dynamics  of  the  same  flow  past  the 
same  body  when  the  latter  is  held  stationary.  The  use  of  a  response 
parameter  given  by  M^/(pLD  ),  as  done  previously  by  others  [31,  without 
(C®  )  could  not  have  correlated  the  response  for  all  smooth  and 
rough  cylinders.  In  fact,  initially,  the  experiments  were  conducted 
with  smooth  cylinders.  It  was  supposed  that  roughness  will  complicate 
the  matters  further.  Contrary  to  this  supposition,  roughness  and  the 
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consideration  of  the  starting  dynamical  conditions  helped  to  delineate 
the  most  important  parameter  governing  the  hydroelastic  oscillations  of 
cyl inders . 

As  noted  in  Chapter  I,  researchers  have  used  different  definitions 
for  the  response  parameter.  One  such  definition  is 

A„  =  -Mr  (8  t r2  S2)  (48) 

r  p  L  02 

This  particular  parameter  has  been  extensively  used  in  correlating 
the  data  obtained  in  steady  flow  about  oscillating  cylinders.  A 
semi-empirical  analysis  by  Iwan  and  Blevins  [1,89]  related  the  relative 
amplitude  in  steady  flow  at  perfect  synchronization  to  A^  by 


- <M4  [Q  3  +  4.52  S.  ,1/2 

4  +  11.94  S‘  4.  +11.94  S* 

r  r 


(49) 


Griffin,  Skop,  and  Ramberg  [61]  proposed  a  least-squares  fit  to 
the  existing  data  as 


-g-  =  1. 29/(1  +  0.43  Ar)3,35 


(50) 


Finally,  Sarpkaya  [10,  73]  theoretically  obtained 


T  ■  9'56  <CL°W”-94,4r2/<cL0)™s  !'/2 


In  which  A  /(C.  °)  is  related  to  R  by 

I  L  III  Id  L/ 


(51) 


(52) 


(cL°) 


rms 


Rp  (8  tt2  S2) 


Thus,  the  variation  of  X/D  with  R  or  A  /(C,  °)  „  in  harmonic  flow 

may  be  compared  with  that  in  steady  flow  through  the  use  of  Eq.  (51). 
Figure  170  shows  such  a  comparison  for  all  the  data  obtained  in  the 
present  investigation.  It  is  rather  surprising  that  not  only  the 
response  in  steady  flow  but  also  that  in  harmonic  flow  is  fairly  well 
represented  by  Eq.  (51).  It  is  of  importance  to  note  that  in  plotting 
the  data  shown  in  Fig.  170  (C^®)^  values  appropriate  to  each  cylinder 
have  been  used  [see  Eq.  (45)]. 

In  summary,  it  has  been  shown  that  the  synchronization  in  harmonic 
flow  occurs  over  a  range  of  Um/fnD  values  from  5.0  to  7.5  (see  Figs. 

47,  51,  55,  59,  63),  with  perfect  synchronization  at  Um/fnD  c  5.6.  The 
amplification  factor  for  the  predominant  harmonic  is  about  2.75,  and 
the  response  is  a  unique  function  of  the  response  parameter  Rp  or 
^r^S-^rms  ^or  smoot^  anc*  rough  cylinders.  The  effect  of  the 
Keulegan-Carpenter  number  and  the  Reynolds  number  on  the  cylinder 
response  is  imbedded  in  the  dependence  of  (C^°)ms  on  the  starting 
dynamical  conditions,  i.e.,  on  the  dependence  of  (C[_°)nns  on  K  and  Re. 
Finally,  the  predominant  harmonic  of  the  response  is  an  order  of  magni¬ 
tude  larger  than  those  of  the  higher  order  harmonics,  which  enable  one 
to  determine  the  response  through  the  use  of  a  relatively  simple  analysis. 
The  modulations  in  amplitude  and  phase  of  the  response  is  quite  well 
predicted  both  by  a  generalized  solution  based  on  the  Fourier  decomposi¬ 
tion  of  the  exciting  force  and  by  a  straightforward  application  of  the 
Duhamel 's  integral. 
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C.  APPLICATION  OF  THE  WAKE-OSCILLATOR  MODEL 

Several  mathematical  models  have  been  proposed  [86-91];  (see  also 
[2,  55,  71,  92]  for  additional  models)  in  an  attempt  to  simulate  and/or 
explain  some  of  the  experimental  results  summarized  in  Chapter  I.  These 
models  do  not  include  the  analysis  of  the  flow  field  and  the  fluid- 
mechanical  arguments  invoked  in  their  evaluation  are  not  altogether 
convincing.  Thus  their  worth  should  be  measured  not  so  much  by  their 
capacity  to  obtain  functional  relations  among  significant  parameters 
that  lead  to  the  basic  understanding  of  the  phenomenon  but  by  their 
ability  to  produce  results  which  are  qualitatively  similar  to  those 
obtained  experimentally. 

The  most  noteworthy  among  the  oscillator  models  is  the  one  proposed 
by  Hartlen  and  Currie  [87]  where  a  van  der  Pol -type  soft  nonlinear 
oscillator  for  the  lift  force  is  coupled  to  the  body  motion  by  a  linear- 
dependence  on  cylinder  velocity.  This  model  is  based  partly  on  a  sugges¬ 
tion  by  Birkhoff  and  Zarantonello  [23],  mentioned  earlier  in  connection 
with  the  estimation  of  the  Strouhal  number,  and  by  Bishop  and  Hassan 
[56]  in  connection  with  their  now  classic  experiments  with  oscillating 
cylinders  in  uniform  flow.  The  model  has  its  roots  in  mechanics  and 
electricity  rather  than  in  the  equations  of  fluid  motion. 

The  pair  of  equations  which  result  from  this  concept  are 

■  aQ  n02  CL  (53) 

K  -  «  ao  h  *  ^L3  +  no2  CL  -  8  V  <54> 
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in  which 


Xr  =  X/D  , 


a 


o 


p  D2  L 

8  7T2  S  2  M 
o 


S  =  f  o/v 

o  o 


T  =  U.nt 


(55) 


Si 


0 


Vf 


n 


where  f  is  the  frequency  of  vortex  shedding  by  a  stationary  cylinder 
in  steady  flow  of  velocity  V  and  SQ  is  the  corresponding  Strouhal  number. 
The  parameters  a  and  y  are  the  van  der  Pol  coefficients  and  B  is  the 
interaction  parameter. 

In  Equation  (54),  first  and  fourth  terms  generate  a  simple  harmonic 
oscillator  of  normalized  frequency  nQ>  the  second  term  (the  so-called 
negative  lift  damping)  provides  the  growth  of  C^,  and  the  third  term 
prevents  unlimited  growth.  The  fifth  term  provides  the  feedback  from 
the  body  motion  to  the  fluid  motion  and  hence  to  the  lift  force. 

Equation  (54)  is  not  the  only  form  which  could  provide  the  desired 
variation  of  C^,  but  it  is  the  simplest. 

Of  the  three  parameters  (a,  y  and  B),  two  must  be  chosen  to  provide 
the  best  fit  to  the  data.  This  is  not  too  unusual  and  all  other  models 
require  some  retrofitting  to  experiment.  The  ratio  a/y  is  related  to 
the  amplitude  of  the  steady-state  oscillation  of  the  lift  force  by  [87] 


CL°  =  *4o73y 


(56) 
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which  is  the  solution  of  Eq.  (54)  for 


Experiments  in  steady  flow  have  shown  that  [56]  in  large-amplitude, 
steady-state,  vortex-induced  oscillation,  the  displacement  and 
exciting  force  have  nearly  sinusoidal  forms  and  oscillate  at  the  same 
frequency  ojc,  close  to  un,  with  a  phase  angle  sj>.  Thus 

xr  =  Xr  sin(fict)  (57) 

and 


CL  =  CL  sin(^ct  +  <t>)  (58) 

Substituting  Eqs.  (57)  and  (58)  in  Eq.  (53)  and  equating  to  zero  the 
resulting  coefficients  of  sin^x)  and  cos(SHct),  Hartlen  and  Currie  [87] 
obtained 


(1  -  n2)  xr  =  aQ  nQ2  CL  cos  0  (59) 

2  ?  X r  =  aQ  nQ2  CL  sin  <p  (60) 

Likewise,  substituting  Eqs.  (57)  and  (58)  in  Eq.  (62),  dropping  terms 

•  3 

arising  from  ,  and  using  Eqs.  (59)  and  (60),  one  has 
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=  nc  n  - 


2a  8?  ft 
o  c 


(l-fic2)2  +  4c2  ft/ 


-1/2 


(61) 


which  shows  that  ft  <  ft  or  f  <  f  . 

co  co 

Hartlen  and  Currie  applied  their  analysis  to  a  comparison  with 
Jones,  et  al . ,  [17],  data  (Re  =  3.6  x  10^  to  1 .9  x  10^,  Xr  =  0.08) 
with  good  qualitative  results.  However,  some  doubt  has  been  raised 
[8]  as  to  whether  the  oscillations  had  reached  synchronization  and 
whether  an  amplitude  of  Xr  =  0.08  was  sufficiently  greater  than  the 
threshhold  amplitude  [44,  62]  necessary  to  bring  about  synchronization. 
Furthermore,  Jones,  et  al,  data  did  not  show  an  increase  in  drag  in  the 
transcritical  range  of  Reynolds  numbers. 

Hartlen  and  Currie  also  used  Ferguson’s  [63]  data  for  the  circular 
cylinder  and  obtained  the  results  shown  in  Fig.  171.  These  results 
are  in  qualitative  agreement  with  those  Shown  in  Fig.  3.  There  are 
however,  some  important  quantitative  differences: 

1.  In  Feng's  data  (Fig.  3)  the  cylinder  is  seen  to  continue  to 

oscillate  at  w  e  u  on  both  sides  of  the  region  outside  the  lock-in 
c  n 

range.  This  feature  is  not  presented  in  Fig.  171. 

2.  The  force  and  response  maximums  occur  in  Fig.  3  at  about  the  same 
Vf-  -alue.  In  Fig.  171,  C^  remains  constant  over  a  broad  range  of 
synchronization,  reaching  ?  maximum  at  a  V  -  value  considerably  smaller 
(relative  to  the  width  of  the  lock-in  range)  than  that  for  X  .  The 
explanation  of  these  discrepancies  and  the  relationship  between  the 
particular  changes  in  the  wake  and/or  the  elastic  system  which  trigger 
the  hysteresis  remained  unresolved. 
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171 


Response  characteristics  in  steady  flow 
as  predicted  by  Hartlen  and  Currie  model  [87]. 
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Fig.  172a  and  Xr  versus  Uf  in  harmonic  flow  as  predicted 
by  Hartlen  and  Currie  model. 
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In  spite  of  its  shortcomings,  Hartlen  and  Currie  model  has  been 


applied  to  the  hydroelastic  oscillations  of  cylinders  in  harmonic 
flow  through  the  use  of  Eqs.  (53),  (54)  and 

CL(t)  =  S_p  si"(V  +  (62) 

The  parameters  a  and  8  have  been  chosen  to  provide  the  best  fit  to  the 
data.  Representative  results  for  a  4  in.  rough  cylinder  are  shown  in 
Figs.  (172a)  and  (172b).  Apparently,  the  judicious  selection  of  the 
two  independent  parameters  provides  a  fairly  good  agreement  between  the 
data  and  the  predictions  of  the  model  near  synchronization.  Calculations 
with  other  cylinders  yielded  similar  results. 

As  noted  earlier,  this  model  does  not  include  the  analysis  of  the 
flow  field  and  gives  no  clues  regarding  the  values  of  a,  y,  or  B. 
Nevertheless,  it  is  worthy  of  consideration  because  of  its  ability  to 
predict  the  general  trend  of  the  data  once  a  and  B  are  chosen  through 
the  use  of  two  data  points  on  the  response  curve.  Additional  attempts 
made  to  improve  the  oscillator  model  by  making  BX^  term  in  Eq.  (54)  a 
non-linear  function  of  Xr  will  not  be  presented  here. 
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VI.  CONCLUSIONS 


The  theoretical  and  experimental  investigation  of  the  hydroelastic 
oscillations  of  smooth  and  rough  circular  cylinders  in  harmonic  flow 
warranted  the  following  conclusions: 

1.  An  elastically-mounted  cylinder  may  undergo  synchronized  oscilla¬ 
tions  when  the  reduced  velocity  Ur  =  Um/fnD  is  in  the  range  of  5.0  to 
7.5. 

2.  Perfect  synchronization,  at  which  the  response  is  maximum,  occurs 

at  U  s  5.6. 
r 

3.  Synchronized  oscillations  occur  at  an  average  Strouhal  number  of 
0.16. 

4.  In  the  region  of  synchronous  oscillations,  the  rms  of  the  lift 
coefficient  and  the  amplitude  of  the  predominant  harmonic  of  the 
normalized  lift  force  are  amplified  by  a  factor  of  about  2.75  relative 
to  that  for  a  stationary  cylinder  in  harmonic  flow  at  the  corresponding 
Keulegan-Carpenter  number  and  Reynolds  number. 

5.  The  Fourier  and  spectral  analyses  of  the  exciting  force  and  response, 
in  the  resonance  region,  show  that  the  exciting  force  is  well  represented 
by  its  three  most  important  harmonics  and  the  response,  by  the  predominant 
harmonic  alone. 

6.  The  analysis  relating  the  governing  parameters  predicts  fairly 
accurately  the  response  of  both  the  smooth  and  rough  cylinders  in  the 
synchronous  region. 

7.  The  measured  instantaneous  values  of  the  cylinder  response  are  in 
excellent  agreement  with  those  predicted  through  the  use  of  the  Duhamel's 
integral . 
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8.  The  wake  oscillator  model  of  Hartlen  and  Currie,  as  applied  to 
harmonic  flow  about  transversely-oscillating  cylinders,  predicts  fairly 
accurately  the  general  behavior  of  the  force  and  response  data  provided 
that  the  two  free  parameters  of  the  model  are  chosen  judiciously. 

9.  Outside  the  synchronization  region  the  cylinders  continue  to 
oscillate,  with  very  small  amplitudes,  at  nearly  their  natural 
frequency. 

10.  The  analyses  and  the  data  presented  herein  are  sufficient  to 
predict  not  only  the  range  of  occurrence  of  hydroelastic  oscillations 
of  elastically-mounted  cylinders  in  harmonic  flow  but  also  to  evaluate 
the  amplitudes  of  the  exciting  force  and  the  response. 
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APPENDIX  A 


PEAK  RELATIVE  AMPLITUDES  AS  A  FUNCTION  OF  K,  Ur,  ff/fnw,  and  fc/fnw 
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APPENDIX  B 


ROOT-MEAN-SQUARE  VALUES  OF  THE  RELATIVE  AMPLITUDES 
AS  A  FUNCTION  OF  K,  Ur,  ff/fnw,  and  yfnw 
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APPENDIX  C 


PEAK  LIFT  COEFFICIENTS  AS  A  FUNCTION  OF  K,  Ur,  ff/fnw,  and  ff/fc 
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APPENDIX  D 


ROOT-MEAN-SQUARE  VALUES  OF  THE  LIFT  COEFFICIENT 
AS  A  FUNCTION  OF  K,  Ur>  ff/fnw>  and  ff/fc 
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versus  K.  for  5  in.  smooth  cylinder,  t  =  0.055. 
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SPECTRAL  ANALYSIS  FOR  A  3  in.  SMOOTH  CYLINDER  FOR  VARIOUS  VALUES  OF  K 
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PHASE  COHERENCE  LOG  SPECTRAL  DENSITY 
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Fig.  135  Spectral  analysis  for  a  3  in.  smooth  cylinder,  c  =  0 
K  ■  60 
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PHASE  COHERENCE  LOG  SPECTRAL  DENSITY 
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PHRSE  COHERENCE  LOG  SPECTRAL  DENSITY 
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SPECTRAL  ANALYSIS  FOR  A  4  in.  ROUGH  CYLINDER  FOR  VARIOUS  VALUES  OF  K 
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PHASE  COHERENCE  LOG  SPECTRAL  DENSITY 


Fig.  144  Spectral  analysis  for  a 
K  =  29 
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a  4  in.  rough  cylinder,  ; 
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PHASE  COHERENCE  LOG  SPECTRAL  DENSITY 


Fig.  150  Spectral  analysis  for  a  4  in.  rough  cylinder,  ;  =  0.030, 
K  ■  43 
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